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"Vol. III, 1952 


Flugbahnen von Leitstrahlraketen mit Gasstrahlsteuerung 
Einfiihrender Bericht 


VR eee) en ee eee 


Von EpuARD PERRET!), ERNST RoTH?), RAYMUND SANGER?) 
und Hans R. VoEttmy!), Ziirich 


1. Einleitung 


ee ee 
¥ - J 


Die Leitstrahlrakete ist dadurch gekennzeichnet, dass sie sich entlang eines 
_ seine Richtung willkiirlich andernden, gebiindelten elektromagnetischen Strah- 
les, dem sogenannten Leitstrahl, bewegt. Wird dabei der Leitstrahl von einem 
_ Radargerat so gesteuert, dass er stets nach einem bestimmten Flugzeug zeigt, 
4 so muss es notwendigerweise zur Begegnung von Rakete und Flugzeug kommen. 
é. Dabei wird angenommen, dass sich die Rakete jederzeit in der Mitte dieses 

Leitstrahls befinde, was einen idealen Servomechanismus der Steuerung voraus- 
setzt. 

Die vorliegende Arbeit befasst sich mit drei verschiedenen Arten der Gas- 
_strahlsteuerung. In den zwei ersten Fallen (1 und 2) wird die zur Steuerung er- 
_forderliche, normal zur Flugbahntangente stehende Kraft durch den gegen die 

_ Flugbahntangente geneigten Gasstrahl erzeugt. Bei der ersten Steuerungsart 
handle es sich um eine iiberstabile Rakete, deren Rumpfachse also stets mit der 
- Bahntangente iibereinstimmt. Der Gasstrahl wird dabei um den Steuerwinkel /, 
iF zur Rumpfachse schraggestellt. Im zweiten Fall sei die Rakete aerodynamisch 

‘indifferent. Der um den Steuerwinkel 8, gegeniiber der Rumpfachse ausge- 

‘lenkte Gasstrahl dreht die ganze Rakete und damit auch den Gasstrahl in die 

- -gewiinschte Schragstellung zur Flugbahntangente. Ein drittes Steuersystem (3) 
: verwendet aerodynamische Tragflachen zur Erzeugung der Querkraft. Diese 

seien fest mit dem Raketenrumpf verbunden, so dass zur Auftriebserzeugung 
die Raketenachse gegen die Flugbahntangente angestellt werden muss. Dem 
um den Steuerwinkel f, ausgelenkten Gasstrahl fallt die Autgabe zu, das riick- 
-stellende aerodynamische Moment zu tiberwinden. 

* Die im Bereich kleiner Anstell- und Steuerwinkel berechtigte Annahme, 
dass zwischen der Drehbewegung der Rakete um Querachsen und der Vor- 
wartsbewegung kein Zusammenhang bestehe, erméglicht die Berechnung der 
~ Flugbahnen, unabhangig von der Wahl des Steuersystems. 


iG 
a 
, 
j 
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An Hand der nun gegebenen Flugbahnen kénnen die Steuerbewegungen der 
erwahnten Systeme berechnet werden. Wirkungstote Raéume ergeben sich dort, 
wo die erforderlichen Steuerbewegungen den Arbeitsbereich der Steuersysteme 
iiberschreiten oder wo andere Grenzen erreicht werden. Die Rakete fallt dann — 
vom Strahl ab und wird zum steuerlosen Flugk6rper. 

Es ware ein zu miihsames Unterfangen, wollten wir das Problem in seiner 
allgemeinsten Form in Rechnung ziehen, und es ist daher angezeigt, von még- 
lichst einfachen Voraussetzungen auszugehen, die aber dennoch erwarten 
lassen, ein zuverlassiges Bild tiber die Bahngestaltung und ihre Bedingungen 
zu gewinnen. In diesem Sinne setzen wir voraus, dass sich das Flugzeug auf 
einer horizontalen Geraden mit konstanter Geschwindigkeit bewege und der 
Abschussort der Rakete senkrecht unter der Bahngeraden liege. In den Rech- 
nungen erscheint neben Flughéhe und Flugzeuggeschwindigkeit der Startzeit- 
punkt hinsichtlich der horizontalen Flugzeugentfernung als Parameter. 

Als Unterlagen fiir die numerischen Rechnungen werden die charakteristi- 
schen Daten eines als Z-Rakete bezeichneten Konstruktionstyps verwendet, 
wobei zur weiteren Vereinfachung der Rechnungen die Abhangigkeit der aero- _ 
dynamischen Koeffizienten von der Mach-Zahl nicht beriicksichtigt wird. 


2. Krafte und Momente, die an der Rakete angreifen 
a) Schubkraft + 


Fiir die durch ein Strahltriebwerk erzeugte Schubkraft gilt 


S=mv,+ 4 (be — pa); (2. 1) 


worin m = dm/dt die pro Zeiteinheit mit der Geschwindigkeit v, (relativ zum 
Diisenausgang) ausgestrahlte Gasmenge bedeutet, q die Diisenaustrittsflache 
und p, bzw. p, den Gasstrahldruck in der Austrittsflache bzw. den atmospha- 
rischen Druck bezeichnet. Es ist zweckmissig, die effektive Gasstrahlgeschwin- 
digkeit vp, einzufithren, in welcher der statische Druckanteil q (be — pq) des 
Schubes miteinbezogen ist, indem wir schreiben 


S=mMv,. (2.2) 


Wir nehmen an, dass die Diisenmiindung drehbar im Raketenende gelagert 
ist, so dass der Gasstrahl bis zu einem maximalen Winkel 


B 1 Bor 


ausgelenkt werden kann. Die Einstellung des Steuerwinkels erfolgt durch das 
Steuergerat der Rakete. 
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b) Gewicht der Rakete 


Mit Riicksicht auf den konstanten Treibstoffverbrauch schreiben wir fiir 
die Masse der Rakete 


t 


M=M,~— | mdt=M,—mt, (2.3) 
0 


wo M, die Anfangsmasse der Rakete bezeichnet. Damit wird das Gewicht der 
Rakete 
G=Mg=Myg-—met. (2.4) 


c) Aerodynamische Krafte und Momente 


Die Resultierende aller aerodynamischen Krafte bei einem auftrieberzeu- 
genden Raketenkérper geht normalerweise nicht durch den Schwerpunkt, 
sondern schneidet die Raketenlangsachse im Abstand / vom Schwerpunkt im 
sogenannten Druckmittelpunkt. 

Die auf den Schwerpunkt bezogene Dyname weist somit folgende Krafte- 
komponenten und Momente auf: 

Den Lujtwiderstand W in Richtung der Bahntangente 


W=-K, ed?v?. (2.5) 


Wegen des ungefahr quadratischen Verhaltens des Koeffizienten Ky als Funktion 
des Anstellwinkels 6 darf fiir kleine 6 der Widerstand W als eine von 0 unab- 
hangige Grésse betrachtet werden. 

Den Auftrieb A senkrecht zur Bahntangente 


A= KK, o870°0= A’ oO. (2.6) 
Das aerodynamische Moment Jp um den Schwerpunkt 
jp] —KpOe 0, (2.7) 
wobei sich Kp aus den aerodynamischen Koeffizienten Ky und K, berechnen 
lasst. Fiir kleine 6 gilt 
Ree (Ry Ken da, Ree hae (2.8) 


Die Grésse / folgt somit leicht aus den beiden Koeffizienten Kr und Ky oder 
auch aus den Momenten und Kraften selbst 


De Krad = Tr (2.9) 


pene A’ 
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Bei aerodynamisch stabilen Geschossen liegt der Druckmittelpunkt hinter 
dem Schwerpunkt, gegensatzlich zu Artilleriegeschossen, die deshalb des Dralls 
bediirfen. Die aerodynamische Stabilitat ist um so besser, je grésser bei gleichem 
Auftrieb die Grosse / ist. 


3. Bewegungsgleichungen 
a) Koordinatensystem 


Da wir, wie eingangs erwahnt, nur geradlinige Flugwege des Zieles voraus- 
setzen, denken wir uns durch den Flugweg und den Leitstrahlsender eine Ebene 
gelegt. Der Leitstrahl bewegt sich bei der Verfolgung des Zieles nur in dieser 
Ebene, und somit tritt eine ideal gesteuerte Rakete nicht aus dieser Ebene 
heraus. Unser Problem der Raketenbahn wird also zweidimensional. 


al 


a 


nl 
r) 
a 


Fig. 3.1 


Koordinatensystem. 


In der eingefiihrten Ebene nehmen wir ein Koordinatensystem an mit dem 
Nullpunkt im Sender. Die x-Achse sei parallel zur Flugrichtung des Zieles, die 
y-Achse senkrecht dazu und nach oben gerichtet (vgl. Figur 3.1). 

Wir kénnten auch ein mit dem Leitstrahl bewegtes Koordinatensystem ein- 
fiihren. Das hatte wohl den Vorteil, dass aus der Rechnung sofort der Abstand 
vom Sender bekannt wiire, aber den Nachteil, dass die Raketenbewegung auf 
ein erdfestes System umgerechnet werden miusste, da einige Koeffizienten Funk- 


tionen der Héhe sind. Ausserdem triaten im bewegten System Tragheitskrafte — 


auf, die eine besondere Rechnung erfordern wiirden. Aus Griinden der ein- 
facheren Darstellung wihlen wir das feste Koordinatensystem. 

Die eingefiihrte (x, y)-Ebene steht senkrecht, wenn das Flugzeug den Sender 
im Zenit tiberfliegt. Diesen speziellen Fall wollen wir unsern Bahnberech- 
nungen zugrunde legen. Die allgemeine Lage wiirde keine wesentlich davon 
abweichende Raketenbahnen ergeben, da allein die in der (x, y)-Ebene liegende 
Gewichtskomponente sich mit der Neigung der Ebene dndert. 
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Wir fiihren folgende Bezeichnungen ein (vgl. Figur 3.1 und O22) % 

A = Neigung des Leitstrahls gegeniiber der x-Achse, 

y = Winkel zwischen Leitstrahl und Bewegungsrichtung, 

gy = Winkel zwischen Bewegungsrichtung und x-Achse, 

0 = Anstellwinkel der Rakete, 

6 = Steuerwinkel, das heisst Auslenkwinkel des Schubvektors gegentiber 
der Raketenachse. 

T= t — alw = Zeit vom Zenit an gerechnet. 


JY Bahn | m eZ 
ell? 
ae ee 


Fig. 3.2 


Krafte, welche an der Rakete angreifen. 


b) Differentialgleichungen der Raketenbahn 


Aus Figur 3.2 ergeben sich in der x- und y-Richtung die Bewegungs- 


gleichungen 
mx= Scos(B+6+ 9)—Wcosp—Asing, (Sik) 


my=Ssin(B+6+ 9)—-Wsing+ 4 cosp—G. (3.2) 
Aus Figur 3.1 ergibt sich ferner die einfache Beziehung 


x= yT. (3.3) 


Die Daten des Flugzeugkurses sind daher in einer einzigen Funktion enthalten 
und gehen in sehr einfacher Weise in das Differentialgleichungssystem ein. 
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Wir bendtigen im folgenden die Differentialquotienten 


=F (yt+yT), (3.4) | 
#=— (2y+yT). (3.5) 


Ferner entnehmen wir der Figur 3.2 die Beziehung 


tgp = =. (3.6) 
Mit den Gleichungen (3. 1) bis (3.6) haben wir alle Beziehungen aufgestellt, 
die zur Berechnung einer Raketenbahn bei idealer Steuerung erforderlich sind; 
sie lassen sich noch etwas vereinfachen. 
Zunachst multiplizieren wir (3. 1) mit cos@ und (3. 2) mit sing und addie- 
ren beide 


mx cosp+my sing = Scos(B + 6) — W — Gsing. (3.7) 


Damit haben wir eine Relation gewonnen, welche den Auftrieb nicht mehr 
enthalt. Entsprechend erhalten wir aus (3.1) und (3.2) mit vertauschter Multi- 
plikation 

—mx sing +mycosp = S sin(B + 6) + A —Geos@. (3.8) 


Durch Einsetzen von (3. 5) in (3. 7) wird 


S cos (B + 0) = 


w es w. i 
m(tep+>T)i+2mey = -—Gtggy. (3.9) 


) cos 
Eliminieren wir noch in (3. 6) x mit (3. 4) 


igp=—- Sia (3.10) 
so haben wir ein Differentialgleichungssystem (3.9) und ek 10) gewonnen, 
welches uns bei gewissen Vernachlassigungen gestattet, die Raketenbahn in ein- 
facher Weise zu berechnen, Diese Vernachlassigungen bestehen in der Annahme, 
f und 6 seien kleine Winkel, so dass cos (6 + 6) ~ 1 gesetzt werden darf. Ausser- 
dem konnen wir bei kleinem 6 den Luftwiderstand W als von 6 unabhangig be- 
trachten. Mit diesen Vernachlissigungen treten in unserem Gleichungssystem 
nur die zwei Unbekannten y und auf; somit lassen sich aus diesen beiden 
Gleichungen allein die Raketenbahnen durch schrittweise Integration ermitteln. 

Fiir die eingangs beschriebenen drei verschiedenen Arten der Steuerung 
belegen wir die Winkel 6 und 6 mit den Indizes 1, 2 und 3. 


ete 
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Die Annahme, £, sei klein bei der, Steuerung mit schraggestellter Brenn- 
kammer, ist nicht ohne weiteres erlaubt, da, wie wir sehen werden, bisweilen 
grosse Steuerwinkel f, des Gasstrahles zur Erzeugung der notwendigen Schub- 
komponente quer zur Flugbahn erforderlich sind. Hingegen ist der Anstellwinkel 
0, sehr klein, weil grosse aerodynamische Stabilitat vorausgesetzt wird. 

Zur Integration muss somit die Gleichung (3. 8) herangezogen werden, welche 
fiir 6, = 0 lautet 

—mx sing + my cosy = SsinB, — Geos. (Sad) 


x kann nach (3. 5) eliminiert werden, so dass die Gleichungen (3. 9), (3. 10) und 
(3.11) ein System fiir die Funktionen y, m und f, bilden. Die Beriicksichtigung 
von f, fiihrt auf eine Verkleinerung der Schubkomponente in Richtung der 
Bahntangente. 

Bei der Steuerung mit schraggestellter indifferenter Rakete kann f, klein 
angenommen werden, sofern der Abstand s (Figur 3.2) zwischen Raketen- 
schwerpunkt und Drehpunkt der Brennkammer geniigend gross gewahlt wird, 
so dass schon eine geringe Auslenkung der Brennkammer die Rakete zu drehen 
vermag. Hingegen kann 6, nicht vernachlassigt werden, da, wie im ersten Fall, 
zur Erzeugung der notwendigen Querkraft eine erhebliche Schragstellung des 
Schubvektors zur Flugbahntangente und daher der Rakete erforderlich ist. 

Zur Integration des Differentialgleichungssystems muss somit auch hier die 
Gleichung (3. 8) herangezogen werden. Wegen f, ~ 0 entsteht wieder die Glei- 
chung (3. 11), mit 6, an Stelle von {,. Die Beriicksichtigung von 6, bedingt hier 
nicht nur eine Verminderung der Schubkomponente in Richtung der Bahntan- 
gente, sondern auch eine Vergrésserung des Luftwiderstandes W, da dieser bei 
merklichen Anstellwinkeln gegensatzlich zu unserer Annahme (2.5) von 6 ab- 
hangig wird. 

Bei der dritten Art der Steuerung sind sowohl , wie 6, zum vornherein klein, 
zum mindesten bei grossen Geschwindigkeiten. Wird namlich die erforderliche 
Querkraft durch grosse Fliigel erzeugt, so geniigen zur Steuerung kleine An- 
- stellwinkel 6. Wird ferner die Rakete sehr wenig stabil angenommen, so reichen 
kleine Steuerwinkel f, des Gasstrahles aus, um der Rakete einen endlichen 
Anstellwinkel zu geben. 

Die dritte Art der Steuerung gestattet somit die beschriebene einfache 


Integration mit (3. 9) und (3. 10). 


c) Erforderliche Querkrafte 


Die normal zur Flugbahn stehenden Krafte nennen wir im folgenden kurz 
Querkrafte. Durch Auslenkung des Gasstrahls entsteht nach Figur 3.2 die Quer- 
kraft S sin(8 + 6). Der Auftrieb als weitere Querkraft entsteht indirekt eben- 


-_ + TT 
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falls durch Auslenkung des Gasstrahls. Die zur Steuerung der Rakete verfiig- 
baren Querkrafte sind somit allgemein 


Q=Ssin(B+o0)+A4. (3.12) 


. 


Aus Gleichung (3. 8) kénnen diese Querkrafte berechnet werden, wenn nach 
der Berechnung der Flugbahn x, y und m als Funktionen der Zeit bekannt sind 


Q=—mxsing+mycosy + Gcosp . (3.13) 


4. Numerische Integration des Differentialgleichungssystems 
der Flugbahnen 


Wir wollen im folgenden vier verschiedene Raketenbahnen wiedergeben, 
denen die untenstehenden Zahlenwerte zugrunde liegen. 


Die Bewegung des Zieles ist charakterisiert durch: 


Flugh6he h = 5km. 
Geschwindigkeit des Flugzeuges w= 200 m/s. 
Horizontaldistanz zur Zeit des Abschusses 
Flug Nr. I II III IV 
a= 10000 m 4000 m 0) — 4000 m. 


In der Rechnung tritt nur das Verhiltnis w/h auf, das heisst, es entsteht die 
gleiche Raketenbahn, wenn das Flugzeug zum Beispiel in doppelter Flughéhe 
mit doppelter Geschwindigkeit fliegt. 


Fir die Z-Rakete nehmen wir folgende Zahlenwerte an: 


Widerstandskoeffizient Ay = 0,086. 

Auftriebskoeffizient Ky = 9,9, 

Raketendurchmesser d = 0,28 m (Aw d? = 0,00674 m2). 
Anfangsmasse my = 25,5 kg* m7! 52 = 250 kg, 

Sekundlicher Triebstoffverbrauch m= 0,204 kg* m=1s = 2 esate 
Schub S = 550 kg* = 5400 kg ms7?. 

Luftdichte auf Starthéhe Op = 0,119 ke*oms 4's*= 117 ko m2. 


In den nachfolgenden Tabellen 1 bis 4 sind die Ergebnisse der numerischen 
Integration zusammengestellt?). 


1) Der Abstand vom Sender ist mit r bezeichnet. Ferner ist zu bemerken, dass die Flugbahnen 
streng unter Annahme der ersten Art der Steuerung mit Beriicksichtigung des Gliedes S cos B (6 = 0) 
in (3.7) berechnet sind. Sie gelten aber mit guter Naherung auch fiir die zweite und dritte Art der 
Steuerung, solange B, kleiner als etwa 25° ist. Dieser Wert wird nur beim ersten Flug tiberschritten, 
wahrend zum Beispiel beim dritten Flug nicht einmal 8° erreicht werden. 


— 
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Tabelle 1 
a =10000m, w = 200 m/s 

t ry ) 

s m m ke* 

0 (0) 0) 230 24° (Oe 

2 Silt 14 234 24°25’ 09302 
4 123 59 240 Zu we 

6 276 136 247 25°40’ 1°40’ 
8 490 250 254 26° 30’ BSS! 
10 762 404 265 272356 soy sys 
12 1090 600 278 29° 4°10/ 
14 1480 842 293 305304 2-207 
16 1910 1130 Sy 322500 6° 40’ 
18 2400 1480 333; 34°40’ 8°10 
20 2930 1880 356 Saf? ty 9°50’ 
22 3490 2330 387 40° 20’ eA® 
24 4090 2840 421 43°50’ 14°25’ 
26 4720 3410 460 48° We AO 
27 5030 3710 485 50°30’ 18°50’ 
28 D350) 4020 509 23P 20°40’ 
29 5670 4340 B55: 57°45’ 22°40/ 
30 5990 4680 625 65° 25° 
Sal 6300 5020 — lit Oe 

Tabelle 2 
a= 4000m, w= 200 m/s 

; | : fs | ° Py yp 

Ss m m kg* 

0 | 0) 10) 158 16°30’ On 

2 26 21 168 17302 il 2i0y 
4 103 87 179 182357 SAOy 
6 230 200 190 19°50’ SOU? 
8 404 365 204 212207 eo 
10 625 580 218 22°45’ O34 
eZ 889 846 230 24° 12°45’ 
14 1195 1162 235 24° 30° oe Om 
16 1540 1520 237, 24°40’ 1 Sra 54 
18 1930 1920 230 24° 20°45’ 
20 2350 2350 218 22°45’ DOSY 
22 2820 2810 200 202504 24°40’ 
24 3330 3280 Neil Pee Zon 25950 
26 3880 3780 154 16° 26° 40° 
28 4500 4280 Wy 1307 272 
30 5170 4800 107 he LO 26°50’ 
32 5910 5330 83 8°40’ 26° 30’ 
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Tabelle 3 
a@=0, w= 200an/s 


ZAMP 


ie 
0 0 0 O GF 
24 24 6 0° 34’ Pai i fs 
96 95 12 Ge 4° 26’ 
219 213 16 1°44’ 622% 
394 376 18 157" SFG" 
624 579 20 ay oh 9° 30’ 
909 819 18 12562 10°39’ 
1253 1094 tS Oe JasSsoe 
1658 1396 11 1°14’ uA ie Fs 
2130 L725 6 0° 34’ 12°40’ 
2660 2070 -— 1 —0° 6’ 127567 
3260 2440 — 7 —0° 45’ UTES lo 
3920 2830 —15 —1°34’ 13°59" 
4650 3230 —22 —2°20’ 132 n" 
5460 3640 —29 —2°58’ 1342 
6 330 4050 —35 —3°35’ 127525 
7280 4480 —40 —4°11’ 12°41’ 
8 290 4910 —45 —4° 38’ 127307 
9 380 5350 —50 —5°12’ 12°16" 
10540 5800 —54 —5° 38’ 1 Ago 
11770 6240 —57 =O 12s, oh Lata) 
13070 6690 —60 —6°20’ us ses bbs 
15900 7600 —65 —6° 52’ Leas 
19020 8500 | —69 —7°16’ 10°33" 
Tabelle 4 
—4000m, w= 200 m/s 


By 


m kg* 

0 0 —158 —16° 30’ 
Lah 80 —136 —14°13’ 
452 305 -118 —12° 28’ 
1054 650 —109 —11° 26’ 
1919 1100 —105 —10° 54’ 
3065 1625 —101 —10°33’ 
4490 2220 —100 —10°22’ 
6190 2860 — 99 —10°17’ 
8170 3540 — 99 —10°17’ 
10420 4250 — 98 —10°14’ 
12940 4980 — 98 —10°11’ 
15720 5730 — 97 -10° 6’ 
18750 6470 — 96 —10° 1’ 
22000 7240 — 95 — 9°48’ 

25600 8000 

29400 8780 


a a ee ee ee ee eee ee S| 


ee on 
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In den Figuren 4.1 bis 4.4 sind die Flugbahnen I bis IV dargestellt. Wie 
wir den Figuren entnehmen, liegt die konvexe Seite der Flugbahn dem Anflug 


des Flugzeuges zugekehrt. Dieser Verlauf ist charakteristisch fiir die Bahneén 
von Leitstrahlraketen. 


w=200M/s 


Fig. 4.1 
Leitstrahlraketenbahn fiir a = 10000 m. 


w=200M/5 


[eee a ee : ieee 
-2000 2000 4000 
Fig. 4.2 
Leitstrahlraketenbahn fitr a = 4000 m. 


‘ 
: 
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40 Jt Os 20s 105 w- 200m/s 


LS 
eS 
Ses 
——~— Xin im 
—8000 -6000 -4000 
Fig. 4.3 
Leitstrahlraketenbahn fiir a = 0. 
60s 0 w-200m/s 
E 
Ass 
a 
-/2000 -10000 -8000 -6000 -4000 | 


Fig. 4.4 


Leitstrahlraketenbahn fiir a = —4000 m. 


5. Berechnung der Steuerwinkel B 
a) Steuerung mit schriggestellter Brennkammer 


Wir betrachten zunachst die erste Art der Steuerung. Da wir in diesem Fall 
grosse aerodynamische Stabilitat voraussetzen, wird 6, = 0, und damit ver- 
schwindet auch der Auftrieb A. Es bleibt somit zur Berechnung von f, aus 
(a, Lz) 


sin, = ©. (5.1) 


— 
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In Figur 5.1 ist der Verlauf von f, fiir die vier berechneten Raketenbahnen 
wiedergegeben. 


° = /reffpunkt 


> Seuerbereich 
3 


rh) 


Big. 5.1 


Steuerwinkel By fiir die vier Raketenbahnen. 


f, erreicht sehr grosse Werte. Wenn wir den Maximalwert £,,. beispielsweise 
20° annehmen, fallt die erste Raketenbahn, fiir welche a = 10000 m, von Anfang 
an aus, da schon beim Start der Grenzwert iiberschritten wird. Auch im zweiten 
Fall wird dieser Grenzwinkel iiberschritten. Hingegen im dritten und vierten 
Fall nicht, so dass hier die Rakete stets dem Leitstrahl zu folgen vermag. 


b) Stewerung mit schrdggestellter Rakete 


Im zweiten Fall geniigen zur Steuerung kleine Auslenkungen der Brenn- 
kammer, weil wir eine indifferente Rakete voraussetzen und nur noch das Trag- 
heitsmoment zu tiberwinden ist. Es handelt sich um eine Differentialsteuerung, 
da die Drehung der Rakete und damit die Anderung der Schubkomponente 
normal zur Bahntangente gesteuert werden muss. Der Winkel f, kann aus der 
Bewegungsgleichung fiir die Drehung um eine Querachse der Rakete berechnet 
werden. Ist © ihr Tragheitsmoment, so ergibt sich aus Figur 3.2 bei aero- 
dynamisch indifferenter Rakete 


O ($+ 6) =—s S fy. (5.2) 


Sehen wir von der erwahnten geringen Anderung der Raketenbahn bei der 
zweiten Steuerung gegentiber der ersten ab, so kénnen wir, weil die Schub- 
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komponente quer zur Bahntangente in beiden Fallen gleich ist, setzen ; 
0» = py . (5.3) 


g und f, sind aus den berechneten Bahnen bekannt. 
Wahlen wir beispielsweise: 
Tragheitsmoment 9 = 20kg* ms? = 196 kg m?. 
Abstand Raketenschwerpunkt — Brennkammerdrehpunkt s = 2m. 
Mit diesen Zahlenwerten erhalten wir fiir die vier berechneten Raketenbahnen 
aus (5. 2) den Verlauf von f, (vgl. Figur 5.2). 


© =Treftpunkt 


Fig. 5. 2 
Steuerwinkel Bs fiir die vier Raketenbahnen. 


Der Steuerwinkel 8, wird auffallend klein (nur einige Winkelsekunden), was 
verstandlich ist, da geringe Auslenkungen der Brennkammer geniigen, um das 
Tragheitsmoment allein zu tiberwinden. Derart kleine Steuerwinkel bedingen 
aber einen empfindlichen Servomechanismus der Steuerung. Der Steuerwinkel 
wird bei dieser indifferenten Rakete unerwiinscht klein, wahrend er bei der : 
iiberstabilen Rakete unerwiinscht gross geworden ist. Es ist somit angezeigt, 
eine gewisse aerodynamische Stabilitat vorauszusetzen, wie dies im dritten Fall 
angenommen wird. 

Handelt es sich jedoch um eine Steuerung im luftleeren Raum, wo prinzipiell 
keine aerodynamische Stabilitat méglich ist, so kommt nur die zweite Art 
der Steuerung in Frage. 


c) Steuerung mit aerodynamischem Auftrieb 


Verwendet man den aerodynamischen Auftrieb von Fliigeln als weiteres 
querkrafterzeugendes Element, so gesellt sich zu den bisher verwendeten 
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Schubgliedern ein geschwindigkeitsabhangiges Glied. In der Startperiode wird 
sich zunachst gegeniiber der Steuerung 2 wenig dndern. Erst mit wachsender 
Geschwindigkeit tibernimmt der Auftrieb immer mehr die tragende Funktion und 
tibertrifft schliesslich bei weitem die bisher wichtig gewesene Schubkomponente. 

Wahrend der Steuerwinkel £, in der Momentengleichung um den Schwer- 
punkt die zentrale Rolle spielt, darf er in den Translationsgleichungen wegen 
seiner Kleinheit neben 6 vernachlassigt werden. Fiir die Querkraft gilt somit 
nach (3.12) 


Os bei a, (5 A’), (5.4) 


welche in der Steuerung 1 durch Q = f, S gegeben war, und durch Vergleich 
erhalten wir die Beziehung 


Bos SrigeAy (5.5) 


Aus der Momentengleichung um den Schwerpunkt folgt der Zusammenhang 
zwischen A und fy (vgl. Figur 3.2) 


A=—B,S =. (5.6) 
Dies ergibt in (5. 5) eingesetzt 
l 
Bs = (63 — fr) ae (5.7) 
Da der Anstellwinkel 6, nach (5. 4) und (5. 5) 


Os = ee Ay Sy (5.8) 
folgt endgiiltig fiir 6, 


1 l 
abel SVESTE Wor 


Zu Beginn der Steuerung ist £; Null, weil A’ als quadratische Funktion von v 
verschwindet. Demzufolge muss der Anstellwinkel beim Start nach (5. 8) gleich 
f, sein. 

Mit zunehmender Geschwindigkeit und damit zunehmendem A’ wird die 
Grésse S/A’ immer kleiner gegenitiber 1, und man erhalt die einfache Beziehung 


y 
b= Pa ts (5. 10) 
Man erkennt nun leicht, dass je nach der Wahl der Gréssen / und s der Steuer- 


winkel f, sehr viel kleiner als 6, gehalten werden kann, wodurch sich der steuer- 
bare Bereich gegentiber-der erstbesprochenen Steuerung wesentlich erhoht. 
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Das negative Zeichen in (5. 10) bringt zum Ausdruck, dass die Brennkam- 
mer gegeniiber der Steuerung 1 in der entgegengesetzten Richtung ausgelenkt 
werden muss! 

In Figur 5.3 sind die Steuerwinkel £6; wiederum fiir die vier Verfolgungsfliige 
dargestellt. Zur Berechnung wurde / = 0,2 m angenommen, was einer sehr kleinen 


2--9llm 


a 4A (zum bernleich) 
© Treftpunkt 


Fig. 5.3 
Steuerwinkel f, fiir die vier Raketenbahnen. 


aerodynamischen Stabilitat entspricht. Mit s = 2m wird //s = 0,1, das heisst, 
die Steuerwinkel f, sind bei grossen Geschwindigkeiten zehnmal kleiner als f, 
[vgl. Formel (5. 10)]. Bei Flugbeginn jedoch ist £3; sehr klein, da seine Funktion 
zunachst nur darin besteht, die Rakete und damit den Schubvektor entgegen 


| Ae Ae 


noch sehr kleinen aerodynamischen Riickstellmomenten in die erforderliche — 


Richtung zu drehen, ahnlich der besprochenen zweiten Steuerung mit indiffe- 
renter Rakete. 

Da fiir alle berechneten Verfolgungsfliige der Steuerwinkel f, nur einige 
Grad betraigt, erkennt man, dass bei dieser Steuerungsart der Wirkungs- 
bereich kaum durch die Begrenzung des Steuerwinkels eingeschrankt werden 
wird. Dagegen ist es nicht einfach, den Abstand / zwischen Schwerpunkt und 
Druckmittelpunkt wahrend des ganzen Fluges klein zu halten, da sowohl der 
Schwerpunkt infolge Treibstoffverbrauches als auch der Druckmittelpunkt 
mit andernder Mach-Zahl wandern. 


d) Wirkungsbereich 


Die hier besprochenen Raketenbahnen stellen wegen der massigen Bahn- 
kriimmungen noch gut erfiillbare Forderungen an die drei besprochenen Steuer- 


—e 
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systeme, mit Ausnahme der ersten Steuerung, die in den Fliigen I und II den 
maximalen Steuerausschlag tiberschreitet. Die zweite Steuerungsart wiirde erst 
dann versagen, wenn die geforderte Querkraft die Gréssenordnung der Schub- 
kraft erreicht, wahrend die auftriebgesteuerte Rakete bei massigen Steuer- 
ausschlagen Querkrafte zu erzeugen vermag, die ein Vielfaches der Schubkraft 
betragen. Die Grenze liegt hier beim héchstzulassigen Auftrieb der Fliigel. 

Die bei der Bahnberechnung nicht beriicksichtigten, aber in Wirklichkeit 
vorkommenden Schwingungen der Raketenbahn um die Leitstrahlmitte und 
der Rakete um ihre Querachsen erfordern zusatzliche Steuerbewegungen und 
Querkrafte, wodurch der Wirkungsbereich weiter eingeschrankt wird. Auch 
besteht die Gefahr, dass der Steuermechanismus nicht mehr richtig auf den Leit- 
strahl reagiert, wenn Flugbahn und Leitstrahl zu grosse Winkel einschliessen. 

Da bei gasstrahlgesteuerten Raketen die Steuerwirkung bei Brennschluss 
erlischt, ist die ausserste Grenze des Wirkungsbereiches durch den erreichten 
Flugweg der Rakete bei Brennschluss gegeben. 


6. Die Rotationsbewegung der Rakete 
a) Drehung um die Querachse 


Die Tatsache, dass die Drehbewegung um die Querachse die Vorwarts- 
bewegung der Rakete nicht wesentlich beeinflusst, gestattete uns, die Raketen- 
flugbahnen unter alleiniger Beriicksichtigung des Widerstandsgesetzes (2. 5) zu 
berechnen. Die Drehung um die Querachse muss erst bei der Berechnung der 
notwendigen Steuerbewegungen beachtet werden. 

Es wiirde den Rahmen dieses Berichtes sprengen, wollte man naher auf die 
Auswirkungen der Drehung um die Querachse eingehen. Bei den aerodynamisch 
stabilen Geschossen sind diese Rotationsbewegungen meist Schwingungen. Diese 
werden durch die formgebundene aerodynamische Dampfung und durch die 
Coriolis-Wirkung des Gasstrahls nur sehr wenig gedampft und bediirfen daher 
der kiinstlichen Dampfung durch die Steuerung. Wegen der Koppelung zwischen 
Drehbewegung um Querachse und Translationsbewegung treten die Dreh- 
schwingungen stets gemeinsam mit Querbewegungen des Raketenschwerpunktes 
auf. Letztere bewirken aber Abweichungen vom Leitstrahl, welche tiber die 
Steuerfunktion die Drehbewegung der Rakete beeinflussen. Derart entsteht 
ein geschlossenes servomechanisches System, zu dessen Untersuchung die Daten 
der hier berechneten «idealen» Flugbahnen die Grundlage bilden. 


b) Drehung um die Langsachse 


Die Drehung um die Langsachse der Rakete spielt fiir die hier berechneten 
Flugbahnen dann keine Rolle, wenn das Geschoss in. bezug auf diese Drehung 
unempfindliche aerodynamische Daten aufweist und wenn die Steuerung fahig 
ist, trotz der Langsrotation die richtigen Funktionen auszuftihren. 


ZAMP III/17 
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7. Bemerkungen zum Servomechanismus der Steuerung 


Es ist in unsern Bewegungsgleichungen angenommen worden, die Rakete 
befinde sich stets in der Mitte des Leitstrahls, was, wie eingangs erwahnt, 
einen idealen Servomechanismus der Steuerung voraussetzt. 

Damit die Rakete bei automatischer Leitstrahlsteuerung (ohne Signale vom 
Boden aus) einen Steuerbefehl erhalt, hat sie sich in einem gewissen Abstand 
von der Leitstrahlmitte zu befinden. Dieser elektrische Steuerbefehl muss so 
beschaffen sein, dass einerseits die Rakete auch bei bewegtem Leitstrahl sich 
mdglichst wenig von der Leitstrahlmitte entfernt und andererseits keine Uber- 
steuerung auftreten kann. 

Befindet sich die Rakete in irgendeinem Zeitpunkt neben der Leitstrahl- 
mitte, so tritt ein Einschwingvorgang ein, bei ruhendem, senkrechtem Leit- 
strahl um seine Mitte, bei schragem oder bewegtem Leitstrahl um eine gewisse 
Gleichgewichtslage. Von der Steuerung wird gefordert, dass dieser Einschwing- 
vorgang moéglichst kurze Zeit dauert. 

Die theoretischen Berechnungen dieser Einschwingvorginge fiihren auf 
Differentialgleichungssysteme hdherer Ordnung mit variablen Koeffizienten, 
die schrittweise integriert werden miissen, wozu programmgesteuerte Rechen- 
maschinen gut angewandt werden kénnen. Auf Grund von Stabilitatsunter- 
suchungen dieser Einschwingvorginge miissen die Steuerkoeffizienten des 
Servomechanismus gewahlt werden, 

Die Schwingungen einer als Massenpunkt betrachteten Rakete um eine 
Gleichgewichtslage beeinflussen unsere berechneten idealen Raketenbahnen 
kaum, da es sich um Bewegungen von nur wenigen Metern handelt gegeniiber 
einer Vorwartsbewegung von einigen Kilometern. 

Wie wir gesehen haben, muss bei der Steuerung durch Auftrieb die Rakete 
wenig stabil sein. Diese kleine aerodynamische Stabilitat wird aber nicht nur 
durch die Wanderung des Druckmittelpunktes infolge variabler Geschwindig- 
keit beeinflusst, sondern auch durch die Wanderung des Schwerpunktes der 
Rakete wegen des Brennstoffverbrauchs. Alle diese Gegebenheiten mitissen in 
der Dimensionierung der Steuerung beriicksichtigt werden. 

Herrn Dr. A. GerBrr, Direktor der Werkzeugmaschinenfabrik Oerlikon, Biihrle & Co., und 
der Contraves AG., Ziirich, méchten wir fiir sein férderndes Interesse an der vorliegenden Arbeit 
herzlich danken. — Fir die Zurverfiigungstellung finanzieller Mittel von seiten der Kriegstechnischen 


Abteilung des Eidgendssischen Militardepartements zur Durchftihrung von numerischen Rech- 
nungen méchten wir unseren verbindlichsten Dank aussprechen. 


Summary 


Assuming a straight coarse of the target, the trajectories of beam rider anti- 
aircraft missiles are computed for different launching intervalls. Three different 
guiding-systems of the jet-guided missiles are discussed by showing the necessary 
stearing movements. 


(Eingegangen: 8. April 1952.) 
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Uber exakte Losungen det Stokes-Navier-Gleichungen inkom- 
pressibler Flissigkeiten bei verinderten Grenzbedingungen 


Von Jakos ACKERET, ETH., Ziirich 


Bekanntlich macht es sehr grosse Schwierigkeiten, Lésungen der hydro- 
dynamischen Gleichungen zu finden, wenn Tragheit und Reibung gleichzeitig 
beriicksichtigt werden sollen. Bei grdsseren Reynolds-Zahlen und schlanken 
Korperformen liefert die Grenzschichttheorie zwar eine praktisch sehr befriedi- 
gende Naherung; bei stumpferen Formen aber setzen Ablésungen ein mit nicht- 
stationdren Wirbeln, die sich einer theoretischen Behandlung fast ganz ent- 
ziehen. 

Nun ist es nur bedingt richtig, wenn man sagt, dass eben die Reibung der 
Gase oder Fliissigkeiten schuld an dieser Schwierigkeit sei. Deren direkte 
Ursache ist namlich das Haften der Flissigkeit an den als starre Kérper be- 
trachteten Oberflachen. Im Falle der Gase weiss man iibrigens heute ziemlich 
gut, dass der Mechanismus der Wechselwirkung eines Gasmolekiils mit einem 
andern («innere» Reibung) sehr verschieden ist von der Wechselwirkung des 
Molekiils mit dem Kristallgitter der Wand («aussere» Reibung). 

Es ist nun mehr diese Grenzbedingung, die die Lésung erschwert, als die 
«cinnere» Reibung. Erst durch das Haften entstehen die grossen Geschwindig- 
keitsgradienten, die zur Grenzschichtbildung fiihren. 

Dass bei anderen Grenzbedingungen unter Umstanden sehr einfache Lésun- 
gen auch bei komplizierten Koérperformen méglich sind, hat schon STOKES 
-erkannt?) 2). Freilich bietet deren praktische Realisierung grosse Schwierig- 
keiten, die wohl nur allmahlich tiberwunden werden kénnen. Trotzdem diirfte 
eine nahere Betrachtung von Interesse sein. Wir wollen uns zunachst mit 
ebenen stationiren Parallelstr6mungen um zylindrische Konturen befassen, 
obwohl die Betrachtungen allgemeinere Giiltigkeit haben. 


1) J. G. Sroxes, Math. Phys. Papers 3, 73. 

2) Nach Fertigstellung der vorliegenden Untersuchung wurde ich in freundlicher Weise auf 
eine sehr interessante Arbeit von Hamer aufmerksam gemacht (Z. angew. Math. Mech. 21, 129-139 
[1941]}), die mir leider entgangen war und die sich auch mit diesem Gegenstand befasst (offenbar 
ohne Kenntnis der Stokesschen Bemerkungen). HAmev hat als Beispiel unter anderem auch den 
Zylinder behandelt [Formel (7) dieser Arbeit]. Es schien mir aber nicht unangemessen, diese Note 
ungekiirzt zum Abdruck zu bringen, da sie, wie ich hoffe, neue Beispiele und Gesichtspunkte 


aufweist. 
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Die Stokes-Navier-Gleichungen : 


ou ou’ 0 (ou av’ op . 
o (use tus) =—5 é + wAu, o (u ix te gp) ee 
reduzieren sich auf die einfachen Euler-Gleichungen ohne Reibung, falls eine 
Potentialstromung vorliegt. 

Dann ist 


Es verschwinden aber auch Au und Av, da zum Beispiel 


O®@ 0 
Ants A(35) =~ (4®)=0. 
Aus den tibrigbleibenden Euler-Gleichungen berechnet sich der Druck genau 
wie im reibungslosen Fall. So ist der Druckwiderstand eines beliebigen Zylinders 
stets Null. Der Auftrieb, falls Zirkulation vorhanden ist, bleibt unverdndert. 

Verandert werden nur die Grenzbedingungen. Zwar werden die Normal- 
komponenten der Geschwindigkeit an der Kérperoberflache immer noch ver- 
schwinden miissen (es sei denn, dass man absaugt). Hingegen miissen tangen- 
tiale Schubspannungen teils mit, teils gegen die Str6mung an der Oberflache 
angebracht werden. Der einfachste Fallist der, dass die Oberflache in sich beweg- 
lich ist und gerade mit der Potential-Oberflachengeschwindigkeit in Richtung 
der Strémung bewegt wird. Eine etwas rohe Realisierung ware méglich durch 
Rollen und Riemen, die natiirlich im allgemeinen von Ort zu Ort verianderliche 
Geschwindigkeiten haben miissen. Falls man dann ideales Haften an dieser 
beweglichen Oberflache voraussetzt, sind Fliissigkeitsbewegung und Grenzbe- 
dingungen der Potentialstrémung vollstandig angepasst. Die schwierige Reali- 
sierung ist wohl der Grund, dass man diese Art Grenzbedingungen kaum 
beachtet, trotzdem der Magnus-Effekt an Zylindern doch deutlich zeigte, dass 
auch mit der dort vorhandenen, sehr unvollkommenen Annaherung schon 
recht beachtliche Resultate erzielt werden konnten}) 2), 

Setzen wir uns aber zunachst tiber solche Bedenken hinweg, so zeigen sich 
sehr einfache Gesetzmissigkeiten, auf die wir hier nun nahe? eingehen wollen. 
Da die Oberfliche bewegt ist und gleichzeitig unter tangentialen Spannungen 
steht, muss Arbeit geleistet werden, trotzdem der Zylinder als Ganzes keine 
Widerstandskraft erfahrt. Die so aufgewendete Leistung ist nun gerade das 
Aquivalent der Dissipation im Fliissigkeitsraume. Die értliche Dissipation ist 


1) Es ist bedauerlich, dass die seinerzeitige iibertriebene Reklame mit dem «Flettner-Rotor» 
das wissenschaftliche Interesse eher verringert hat. 

2) Im Jahre 1918 schlug RAytercu in einer kurzen Note (Sci. Papers 6, 552/553) vor, rotierende 
Zylinder bei zweidimensionalen Diffusoren zu verwenden. 


2 a 9 ee teen 
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bekanntlich im ebenen Fall: 


ey ue (eS) (= )' ie Ou \ 2) 
y oe Ox es Oy & oe oa le 
Uber das ganze Stroémungsfeld integriert, ergibt sich bei wirbelfreier Stromung 
die Gesamtdissipation D (mit du/dx = — 0v/d0y und du/dy = 0v/Ox) zu 


D=bf | daxdy=pof [ {4 (se) +4 (Se) ax ay 
=a] fle yt 4 (Roca 
Schreibt man 


p- ans f(sey+ (I+ (yy ee. 


so ergibt eine einfache Anwendung des Greenschen Satzes (da ja Au und Av 
verschwinden) 
ae 


iiber alle Grenzflachen integriert. Das ist auch gleich (mit c? = u? 4+ v?) 


Da=-pb fo ds=— 2 fo, Fe ds. (1) 


Nun ist aber die Schubspannung in der reibenden Fliissigkeit 


Ou Ov 
rau (se + sy) 
und an den Oberflachen: 
OC, 
le ty, 


in Richtung der Bewegung wirkend. 
Wir erhalten somit fiir die aufzuwendende Leistung L: 


D=L=—b | t0,ds (2) 


ein Ergebnis, das auch anschaulich erkannt werden kann. 
Fiir die praktische Berechnung lasst sich dieser Ausdruck noch weiter ver- 


einfachen : 
Bei Potentialstro6mung ist namlich: 


0c, _ ny 


an as." 
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Beachtet man, dass (Figur 1) 
dc,=c¢c,d0 und Rdd=—ds, 
so folgt in der Grenze 
\ a es 
ip Mae 8% 
also 
2 uc, 
To = — R ; 
Fig. 1 
Hodograph. 
so dass die aufzuwendende Leistung sich schreibt: 
oa 
L=2ub| Fads, (3) 


das Integral tiber den Umriss der Kontur erstreckt. Der Anteil der unendlich 
weit entfernten iibrigen Grenzen ist verschwindend. 


Man kann einen aquivalenten Widerstand W* formell einfiihren durch: 
9 RAR 
L=W*V, We = 28? [Sag 
oder dimensionslos mit einer charakteristischen Lange L 
we mn c2 
Ki x = : ead 
eS ga PEPE Oe pT E Ty es 


und mit Einfiihrung der Reynolds-Zahl Re — oV Llu 


Ae PENS? ae + 
* i ar Riga 
faite call (+) FoR (4) 


ee 


ZAMP ~ 
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c% geht daher im ebenen Fall umgekehrt proportional zur Reynolds-Zahl, 
nimmt also sehr schnell ab. 
Da fiir 4hnliche Konturen das Integral J gleich gross ist, folgt, dass die 
Leistung 
= g tee 7 
bec Vb = 22d J 


nur von der Geschwindigkeit, nicht aber von der Dimension der Kontur ab- 
hangt. Dieses zunachst etwas merkwiirdig erscheinende Ergebnis erklart sich 
anschaulich daraus, dass die Geschwindigkeitsgradienten und damit die Schub- 
spannungen umgekehrt proportional zum linearen MaBstab der Kontur gehen. 
Formel (4) ist fiir alle Reynolds-Zahlen giiltig. 

Da der Kriimmungsradius FR vielfach eine recht komplizierte Funktion der 
Bogenlange ist, ist unter Umstanden eine dritte Schreibweise von Vorteil, die 
die Integration in die Hodographenebene verlegt (Figur 1). 

Es ist dem Betrage nach P 
cs ds 2 ; 

See dd ; 
das ist aber 2d’, wo dF das schraffierte Flachenelement im Hodographen 
bedeutet. 

Somit stellt 


L, »=4yb [dF =4ydF-, (5) 
1 
die Leistung dar, die auf das Oberflachenstiick 1, 2 entfallt. 


Beispiele 
1. Kreiszylinder im Parallelstrom mit Zirkulation (Figur 2) 


Hier ist, da R=const, Formel (3) 
geeignet. Es ist c, = 2 V (sing + sing*), 
wo y* der Winkel ist, der zu den Stau- 
punkten gehért. Er hangt mit dem Auf- 
triebsbeiwert c, (bezogen auf Zylinder- 
durchmesser) wie folgt zusammen: 

sin p* = * ‘ 
Formel (3) ergibt nach elementarer 
Rechnung: 


c2 } Fig..2 
72 ag y 
L= 8x fad tae {1 a 8 72 J (6) Zylinder im Parallelstrom mit Zirkulation. 
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beziehungsweise 
16 x Be 
* a 4 & 
Cw Rea {1 + eal J? (6a) 
wo 
V2R 
Reg= a, 
ewe Lu 


Fur symmetrische Anstrémung (keine Zirkulation, also auch c, = 0), folgt: 
CL=8n avo: (7) 


Bei sehr grosser Zirkulation wachst die Leistung rasch, bei p* = 90°, c, = 42 
ist sie beispielsweise dreimal so gross wie im symmetrischen Fall. 


2. Schwach gewolbtes Wellblech (Figur 3) 


Fig. 3 
Uberstrémtes Wellblech. 


Seas h 
Yo = hsin— —, wo —<l. 
Das Potential ist 


5 22 
&=Vx—Vheos=** ee 
die Geschwindigkeit an der Oberfliche: 
Re a ees x 
ee V {1 + ; "sin “FAI, 
die Kriimmung: 
1. Aath Bec x 4 2? 
R= 2 ewes ed) 


Nach Formel (3) wird fiir eine volle Periode 


L=2yb| Sdx=16 0% yV2o (8) 
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Der aquivalente Widerstand ist 


3. Elliptischer Zylinder (ohne Zirkulation) (Figur 4) 


Hier ist es vorteilhaft, nach For- A 
mel (5) zu rechnen, das heisst zuerst 
den Hodographen zu bestimmen. Wir 
bentitzen die konforme Abbildung 
des Kreises auf die Ellipse. 

Die Abbildungsfunktion ist 


S, 1 a* 
biaelerh Sa) 


Fir die Halbachsen erhalten wir uy ., 
; v 
ee Ce 


2 
a=F(1 +3), 


ee 


v2 
Y a? Fig. 4 
B Fis (1 =i vr) Konforme Abbildung des Kreises auf die Ellipse. 
Das Verhaltnis o = B/a ist also = (1 — a?/r?)/(1 + a?/r?), damit 
a 1—o a? 2 
i Se Tae ek ip eae * 


der Radius des abzubildenden Kreises ist damit gegeben. 
Die Geschwindigkeit im Aquatorpunkt A ist vor der Abbildung 2 Vz, = V, 
nach der Abbildung 


- V V 
Mr ah some oe ce i We a? \ | 
| dz |4 call =i) 
Daas, = 27; folet 
Zave 
1 = gr = V+) 
ae 


Man kann nun leicht zeigen, dass der Hodograph fiir ein Ellipsenviertel der 
Halbkreis mit dem Durchmesser ¢, | ist. 


Allgemein ist namlich: 
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und mit z=re'” 


iy as ae 
eee 4\1 oa = Ans cos 2 7}, 
a” sin2¢ re 
Cys eee 
PCOS a1 


Die Geschwindigkeit c, ist gegeben durch 


ee ee 2Vsng _ y 
dai otra Waist Perna 
|/: = < irs cos 2 p 


Berechnen wir: cos} = sin (p+ y) , so finden wir nach einfacher Rechnung: 
1 z a? 
cos & = 2r* sim ~ (1 + ar)» 
so dass also die Beziehung besteht: 
cs= cos? V (1+ o) = Csman COSY - 


Damit ist gezeigt, dass der Hodograph fiir jede umstrémte Ellipse ein Kreis ist. 
Die erforderliche Leistung kann nach Gleichung (5) berechnet werden: 
Fir die Viertelellipse ist 


ZL =4 pb Eyamy, oder L=22pnV2b (1 +0)? (9) 
Ist 6 = 1 (Kreiszylinder), so folgt 
L=8xuV?d, 


in Ubereinstimmung mit Gleichung (7). Interessant ist, dass fiir o = 0 (ebene 
Platte parallel angestrémt) ein endlicher Wert herauskommt: 


L priate = 2% VPD, (9a) 


Der Grund dafiir liegt in der sehr grossen Dissipation in der Nahe der Stau- 
punkte. Die Lange der gekriimmten Linienelemente nimmt zwar ab ; dafiir 
wird der Kriimmungsradius kleiner. 


Bei Querstellung der Platte ist ¢ = oo, und damit wird die theoretische 
Leistung auch unendlich gross, 


—— 


- uf i ° 
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4. Andere Profilformen 


Man kann an Hand des Hodographen sehen, dass konvexe Ecken die Ar- 
beitsleistung unendlich, konkave aber einen endlichen Wert dafiir ergeben 
(Figur 5 B und 5 A). 


Fig. +) 


Hodographen bei ein- und ausspringenden Ecken. 


Sind in a und c irgendwelche endlichen Werte von c, vorhanden (Figur 5 4), 
so ist die Hodographenflache endlich — im Falle der Figur 5 B aber unendlich. 
Von unserem Standpunkt aus muss man also konvexe Ecken vermeiden. 
Fiir eine Sichelform kann man den Hodographen sofort schatzungsweise 
angeben. Man erkennt, dass die Hodographenflache bei diinner werdender 
Sichel gegen Null geht (im Gegensatz zur unendlich diinnen Ellipse) (Figur 6). 


6 =) C 


Fig. 6 
Hodograph des Zweieckprofils. 


Es ist nicht verwunderlich, dass derartige Unterschiede vorkommen, da ja 
hier der kriimmungsmildernde Einfluss der Grenzschichten weefallt und die 
Dissipation sehr hohe Werte erreichen kann. 
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Ein besonders einfacher Fall entsteht, wenn Begrenzungen mit konstanter 
Geschwindigkeit vorkommen (die auch leichter zu realisieren waren). 

Denken wir uns die gestrichelte Freistrahlgrenze in dem von RIABOUCHINSKY 
und DEMTSCHENKO!) behandelten Fall (Figur 7) materiell, aber auch in sich 
beweglich, so kann die Leistung so- 
fort angegeben werden, wenn c,/V 
aus der konformen Abbildung gefun- 
den ist. Der Hodograph fiir einen 
Quadranten ist in Figur 7 unten skiz- 
ziert. 

Die Leistung ist nach Formel (5) 


J a4ubFynca=4uo mcf 
10 
L=4nucb. ( ) 


Man kann die Leistung in zwei 
Schritten auch nach der Formel (3) 
rechnen: 

Von b bis c ist R = ov, der Anteil © 
an L also Null. 

Vonc bis dist c, = cy = const, also: 


Fig. 7 
J d 
Hodograph einer geschlossenen Freistrahlstré- ds 
mung nach RraBoucHinsky und DEMTSCHENKO. = 2 lt b CF f- rz 


c 


Da aber ds/R = dd und / dd = 2/2, so folgt wiederum (10). 


5. Kugel 


Wir haben bisher ebene Probleme behandelt. Es ist nicht schwer, auch : 
andere — etwa rotationssymmetrische Formen — zu untersuchen, zum Beispiel 
die Kugel. 

Fiir den Meridianschnitt gilt wie friiher: 


T=—2y a 
Mit 


oped 
= 7 Vsing und df=22 Rsing Rdp 


1) RIABOUCHINSKY und DremrscHENKO, Verh, 3. int. Kongr. techn. Mech., Stockholm 1930, 
S. 159-165. 
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ergibt sich 


L=— | %e,df=12npV?R. (11) 
Definieren wir 
L V2R 
a : = : AL Ve 
Cy = (0/2) Ven B® und. Ke; = cae 
so ergibt sich 
48 
OF Tes (11a) 


Vergleich mit normalen Grenzbedingungen 


Ziehen wir etwa die Kugel heran, so zeigt unsere Formel eine sehr nahe 
Verwandtschaft mit der fiir langsame Bewegungen (Re, < 1) giiltigen Stokes- 
schen Formel. Der Betrag des aquiva- 
lenten Widerstandes ist genau um den 
Faktor 2 grésser als der Stokessche. % 
Dass dieselbe Abhangigkeit von Re vor- 
handen ist, kann man sich durch die Be- 
trachtung der Geschwindigkeitsgradi- 
enten klarmachen, die in beiden Fallen 
von der Gréssenordnung V/Rsind.Na- 0" 
tiirlich hatte es unter diesen Umstanden 
keinen Nutzen, bewegliche Oberflachen 
zuverwenden.(Beispielsweise wiirdebei 7 
Bakterien, wo mitbewegte Oberflachen 6, = a. 
vielleicht noch denkbar waren, damit si ee : 

; ; : j okes) 
kein Gewinn entstehen.) Die logarith- — g-|__ sae 
mische Auftragung, Figur 8, zeigt sehr 
deutlich, dass der Nutzen der vollen 
Mitbewegung bei grossen Reynolds- __, 
Zahlen liegt. Beispielsweise wiirde der fe 
aquivalente Widerstandswert bei Fe 
= 480000 (Gréssenordnung der kriti- es 
schen Reynolds-Zahl ftir Umschlag J wo’ w°? 10 w? 0? 0° 1° 10°Re® 
der Grenzschicht von laminarer in tur- 
bulente Str6mung) auf 1/10000 sinken, boas ; 

is : cs : : Kugelwiderstand bei normalen Grenzbedingun- 
wahrend die gewohnlichen Grenzbedin- gen und aquivalenter Widerstand bei bewegter 
gungen ~ 1/5 liefern. Ahnliche Verbes- Oberflache. 
serungen findet man auch im zwei- 
dimensionalen Fall, wo praktisch Reynolds-Zahlen von Millionen auftreten, 
so dass die Leistung, die zur Uberwindung der Reibung erforderlich ist, 
praktisch verschwindet. 


Kugel- 
Messungen 
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Stabilitat 


Angesichts dieser Méglichkeiten drangt sich, ganz abgesehen von techni- 
schen Schwierigkeiten, die Frage auf: Ist die hier vorausgesetzte laminare 
Strémung bei mitbewegter Oberflache auch bei grossen Reynolds-Zahlen stabil ? 

Betrachten wir zundchst ein einfaches Beispiel: die Kanalstromung mit 
fester bzw. mitbewegter Wand (Figuren 9a, b, c). Im zweiten Falle kann man 


| [| 


=> 


—> 
teste Wand bewegliche Wand Stérung im 
Relativsystem 
Fig. 9 


Zur Stabilitat der Strémung bei mitbewegter Wand. 


natiirlich sogleich auf Ruhe transformieren, indem man in ein mitbewegtes 
Koordinatensystem geht. Nun gibt aber eine beliebige, an der Wand ver- 
schwindende kleine Stérung (Figur 9c) in sonst ruhender Flissigkeit bei 
ruhenden Wanden keinerlei Anlass zu Anfachungen. Es wiirde ja eine Energie- 
quelle fehlen, so dass mindestens in diesem Falle die Strémung absolut stabil 
ist. Schwieriger diirfte die Frage im allgemeinen Fall zu beantworten sein. Doch 
kann man etwa folgendermassen argumentieren. Bei gebrauchlichen Anwen- 
dungen, wie etwa Flugzeugen oder grossen Objekten im Windkanal, auch bei 
Wasserstr6mungen im grossen Massstab, hat man noch nie Turbulenz im freien 
Flissigkeitsraum, wo Potentialstr6mung herrscht, entstehen sehen. Stets kam 
die Turbulenz von den Wanden. Da wir Haften relativ zur bewegten Wand 
angenommen haben, wird sich eine beliebige Stérung in der wirklich freien 
Str6mung und bei mitbewegter Wand etwas verschieden verhalten miissen 
insofern, als dort die Strémung relativ zur Wand verschwinden muss. Es kén- 
nen sich also dort auch Grenzschichten ausbilden. Es wire ein Beweis notig, 
dass diese nicht zu Turbulenz fiihren kénnen. Die bisherigen Versuche zum 


Magnus-Effekt sprechen eher zugunsten der Laminarhaltung der Aussen- 
stromung. 


— ow 


ag 
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Teilweise Mitbewesung 


Die genaue Realisierung der Grenzbedingungen (mitlaufende Wande) wird 
im allgemeinen so schwer sein, dass mit einer idealen Erfiillung kaum zu rech- 
nen ist. Um so mehr k6énnte eine teilweise, eventuell auch stiickweise Mitbe- 
wegung von Bedeutung sein. Freilich treten dann Grenzschichten auf, und die 
schone Einfachheit unserer bisherigen Uberlegungen verschwindet. Es diirfte 
sich aber meines Erachtens lohnen, die Methoden der Grenzschichttheorie auf 
teilweise mitbewegte Wande auszudehnen. Schon mit einfachen Abschatzungen 
kann man erkennen, dass zum Beispiel bei Druckanstieg die Ablosung weit- 
gehend vermieden werden kann, wenn die Wand auch nur einen Bruchteil der 
potentiellen Oberflachengeschwindigkeit hat. Schwieriger diirfte die Frage der 
Stabilitat sein, da man mit stark gekriimmten Geschwindigkeitsprofilen rech- 
nen muss. 


Summary 


As STOKES pointed out long ago there exist exact solutions of the hydro- 
dynamic equations if the boundary conditions are modified. If it can be arranged 
that the «solid» boundaries are actually moving with the fluid velocity then really 
simple solutions are obtained. Several examples are worked out. As turbulence is 
unlikely to occur under these conditions the power to move bodies in a viscous 
fluid at high Reynolds’ numbers would become extremely low. It seems that even 
a partial fulfilment of these changed conditions would be of practical interest. 


(Eingegangen: 13. Dezember 1951.) 


Uber die Behandlung von Mehrleitersystemen mit transversal 
elektromagnetischen Wellen bei hohen Frequenzen 


Von Hans J. VON BAEYER und RONALD KNECHTLI, Baden!) 


1. Einfihrung 


Die Fortpflanzung elektromagnetischer Wellen langs eines Systems parallel 
verlaufender Leiter spielt in der Technik der kurzen und ultrakurzen Wellen 
eine grosse Rolle. Wird die Welle dabei innerhalb eines von einem einzigen Lei- 
ter berandeten Gebietes gefiihrt, so handelt es sich um eine «Hohlleiter»-Welle, 
deren dampfungsarme Ausbreitung bei gegebenen Leiterabmessungen erst ober- 
halb einer bestimmten Minimalfrequenz (der Grenzfrequenz) méglich ist und 


1) NG Brown, Boveri & Cie., zur Zeit Montreal (Kanada) bzw. Massachusetts Institute of 
Technology, Boston (USA.). 
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bei welcher stets die Langskomponente entweder des elektrischen Feldes oder 
des magnetischen Feldes von Null verschieden ist. Im ersteren Fall spricht man 
von einer transversal magnetischen (TM.), im letzteren Fall, von einer trans- 
versal elektrischen (TE.) Welle. Sind dagegen mehrere Leiter an der Wellen- 
ausbreitung beteiligt, so ist neben den TE.- und TM.-Wellen noch eine TEM.- 
Welle méglich, bei welcher die Langskomponenten sowohl des elektrischen als 
auch des magnetischen Feldes iiberall Null sind. Da bei dieser Welle keine 
Grenzfrequenz existiert, bleibt sie bei tiefen Frequenzen, bei denen weder TE.- 
noch TM.-Wellen auftreten kénnen, als einziger und bedeutungsvoller Wellen- 
typ bestehen. 

Wahrend nun die Eigenschaften der TEM.-Wellen auf Koaxial- und Par- 
alleldrahtleitungen allgemein bekannt sind, wurden Mehrleitersysteme mit 
unsymmetrischen Anregungs- und Abschlussbedingungen sowie mit beliebigen 
Querschnittsverhaltnissen der einzelnen Leiter in der Kurzwellentechnik bisher 
wenig behandelt. In der Kabel- und Freileitungstechnik ist die Theorie solcher 
Mehrleitersysteme schon lange bekannt und Gegenstand zahlreicher Unter- 
suchungen [1]#), jedoch ist deren Anwendung auf Kurzwellenprobleme mit den 
speziellen, dort angebrachten Vereinfachungen und Randbedingungen noch 
wenig untersucht worden. Da sich nun die Falle mehren, in denen eine strenge 
Behandlung des Mehrleiterproblems in der Kurzwellentechnik notwendig wird, 
und da die Beherrschung solcher Falle neue technische Méglichkeiten erdffnen 
wird, werden im folgenden, von den Maxwellschen Gleichungen ausgehend, die 
verallgemeinerten Leitungsgleichungen in einer zur Anwendung in der Kurz- 


wellentechnik geeigneten Form aufgestellt und zwei Anwendungsbeispiele er- — 


lautert. Das eine Beispiel betrifft einen Richtungskoppler, das andere die 
gefaltete Dipolantenne, die in erster Naherung als Dreileitersystem aufgefasst 
werden kann. 


2. Allgemeine Gleichungen eines n-Leiter-Systems 


2.1. Feldgleichungen 


Als Mehrleitersystem bezeichnen wir eine Anzahl parallel verlaufender Lei- 
ter beliebigen Querschnitts. Die Summe der Einzelleiterladungen pro Langen- 
einheit sei in jedem Querschnitt gleich Null, was erfiillt ist, wenn einer der 
Leiter die anderen in Form einer auf Potential Null liegenden Abschirmungs- 
hiille umgibt oder das System iiber einer unendlich ausgedehnten leitenden 
Ebene (Erde) verlauft. Die Langsachse sei mit z, die Koordinaten der Quer- 
schnittsebene mit x und y bezeichnet (Figur 1). 


1) Die Ziffern in eckigen Klanimern beziehen sich auf das Literaturverzeichnis auf Seite 285, 
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Setzen wir nun die Existenz eines Systems von TEM.-Wellen voraus, so 
muss tiberall ; 


Ee oF SH S0 (1) 


Fig. 1 


Allgemeines Vierleitersystem, in welchem einer der Leiter als Abschirmung fur die anderen dient. 


sein. Somit werden die Maxwellschen Gleichungen: 


OE OH, 
ier aire (2) 
LIA eiyop OH, 

rot £ = ot a ha an , (3) 
rot, E =0, (4) 

OH, 0E 
He ok, 40 te, 8 

OE 

rotH =ok +e-5 oy = E, Bee (6) 
rot,H =0. (7) 


Aus (4) und (7) sehen wir, dass E und Fi in der (x, y)-Ebene von Poten- 
tialen ableitbar sind, so dass man schreiben kann: 
i = —gradp, H = —grady. (8) 
Aus den vier anderen Gleichungen erhalt man bei Beschrankung auf die har- 
monische Lésung und fiir den ausserhalb des Leitermediums liegenden Raum, 
in welchem wir der Einfachheit halber o = 0 annehmen, 


eb 1 oF EL se (alee OEE (0) 
Oz2 v2 «Ot? ” Oz? ov Ot? 


ZAMP III/18 
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mit v? = 1/(e w). Dies ist zusammen mit Gleichung (1) der Ausdruck fiir unge- 
diampfte TEM.-Wellen, deren Fortpflanzungsgeschwindigkeit unabhangig von 
der Konfiguration der Leiter gleich der Lichtgeschwindigkeit ist. 

Physikalisch erfiillt jedes verlustlose Mehrleitersystem mit parallel verlau- 
fenden Leitern beliebigen Querschnitts die Bedingungen fiir das Auftreten 
solcher Wellen, jedoch ist die TEM.-Lésung nur dann die einzige Lésung der 
Maxwellschen Gleichungen, wenn die Leiterabsténde und Querschnittsdimen- 
sionen klein gegen die Betriebswellenlange sind. Diese Einschrankung ist fiir 
Systeme mit nur zwei Leitern wohlbekannt und sei hier ohne strengere Begriin- 
dung vorausgesetzt. 


2.2. Sirime und Spannungen im Mehrleitersystem 


Da die Feldgréssen E und H nach (8) von skalaren Potentialen ableitbar 
sind, gelten in der (x, y)-Ebene die statischen Feldgleichungen. So besteht 
zwischen den Ladungen q (pro Langeneinheit) der Einzelleiter und deren Poten- 
tialen # der bekannte [2] lineare Zusammenhang 


91 = Cy Py + Cig Pg t * <>, 


(Cin = Cri) 
92 = Coy Py + Cop fat ---, 
der sich in Matrizenform schreiben lisst : 
q=cp. (10) 
Hierin sind die ¢,, cy,, ... als die elektrostatischen Induktionskoeffizienten 


definiert, welche sich aus der geometrischen Konfiguration des Leitersystems 
ergeben und im vorliegenden Fall entsprechend den Ladungen auf die Langen- 
einheit des Systems bezogen sind. Fiir den praktischen Gebrauch eignen sich 
als unmittelbare Messgréssen besser die sogenannten Direktkapazitaten zwi- 
schen den einzelnen Leitern (wieder pro Langeneinheit des Systems), welche 
aus (10) abgeleitet werden kénnen [2] und mit den Induktionskoeffizienten 
folgenden Zusammenhang haben: 


— 0:3 = Cy, ftir “ee Re (3,4 = 1, 2;3,.4;, 9) 


Dabei ist die Erde oder der das System abschirmende Leiter mit dem Index 0 
bezeichnet und sein Potential gleich Null gesetzt. Dieser Nulleiter, der im rein 
elektrostatischen Fall auch im Unendlichen liegen darf, ist bei den uns inter- 


<= mene, gmt 
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essierenden Hochfrequenzanwendungen wichtig, damit das ganze System inner- 
halb eines gegentiber der Wellenlange kleinen Querschnittgebietes nach aussen 
hin abgeschlossen ist. 


Weiterhin ergibt sich ein Zusammenhang zwischen den Potentialen # der 
Einzelleiter und den auf diesen fliessenden Strémen J iiber den allgemein 
giiltigen Kontinuitatssatz. Dieser verkniipft zunachst die zeitliche Ladungs- 
anderung des 7-ten Leiters mit der Stromanderung langs dieses Leiters nach der 


Beziehung: 2 : 
ag t= Sion Ye ey 


Setzt man dies in (10) ein, so erhalt man ein neues Gleichungssystem, das 
wieder in Matrizenform geschrieben werden kann: 


0 0 
ie geo ue) 


und welches die Abhingigkeit des Stromes auf jedem Einzelleiter von der 
Potentialanderung saimtlicher im Leitungssystem vorhandenen Leiter angibt. 

Die Existenz ungedimpfter TEM.-Wellen im Raume ausserhalb des Leiter- 
mediums nach den Gleichungen (9) hat zur Folge, dass entsprechend den Feld- 
grossen E und # auch die Potentiale und Stréme aller Einzelleiter separaten 
Wellengleichungen gehorchen. Integriert man naémlich gemass (8) die elektri- 
sche Feldstirke von der Oberflache sy des Nulleiters bis zur Oberflache s; des 
i-ten Leiters und beachtet, dass das Potential p, des Nulleiters gleich Null ist, : 
so erhalt man das Potential des i-ten Leiters 


p= — | Bas. 


Fiihrt man diese Integration auf beiden Seiten der Wellengleichung (9) durch, 
so entsteht die Wellengleichung 


0? i ab 0? 
Oz2 p= v2 Of Pi 


mit der allgemeinen, periodischen Loésung 
pret eB eka ts), (Isai 


worin y = w/v = 2n/A gesetzt ist. ! 

Auf aibnliche Weise fiihrt die Verkniipfung des Stromes J; mit dem Magnet- 
feld H an der Oberflache des i-ten Leiters nach Einsetzen in (9) zu einer Wellen- 
gleichung a 
gat =o on 
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mit der allgemein periodischen Lésung 
Je=Dieil? 7". 6 eg eae. (13b) 


Die Konstanten A, B, D und F stellen dabei die Amplituden der hin- und riick- 
laufenden Potential- und Stromwellen auf dem i-ten Leiter dar. Die Ampli- 
tuden D und F lassen sich nun unter Benutzung von (13a) und der Verkniip- 
fungsgleichung (12) auf die Amplituden A und B zuriickfiihren, wobei man 
wieder vorteilhaft von der Matrizenschreibweise Gebrauch macht, da nicht 
mehr nur der 7-te Leiter, sondern auch alle iibrigen Leiter mitberiicksichtigt 
werden miissen. Fiihrt man statt der Potentiale die Spannungen U ein, 
welche den Potentialdifferenzen gegeniiber dem Nulleiter entsprechen und 
setzt zur Abkiirzung die Spannungsamplituden 4 e/®! = U’ und Bei®t = Us 
so erhalt man schliesslich das Gleichungssystem: 


1 4-12 NM +7yz 
U=U'e 1v + U"etl* 


(15) 
J=we {Oe Ulett 
welches die grundsatzliche Spannungs- und Stromverteilung auf dem Mehr- 
leitersystem als Uberlagerung von hin- und riicklaufenden Spannungswellen 
angibt. 
In vielen Fallen liegen als Randbedingung die Stréme und Spannungen zu 
Anfang oder Ende des Leitungssystems vor. Dann ist es vorteilhaft, den Glei- 
‘chungen (15) eine andere Form zu geben, in welcher statt der Konstanten U’ 
und U" die Spannungen U, und die Stréme J, am Leitungsanfang (z = 0) 
auftreten. Aus (15) ergibt sich namlich 


U,=U'+U", 
Jo = ve{U'— u"\ 
und hieraus: 
PRES Finis Se 
U"= = |U,— = et Jo. 


Setzt man dies wieder in (15) ein, so erhalt man: 


U= U, cosy z — j—etJysinyz, 
(15a) 
J=Jycosyz—jueUjsinyz. 


Das Gleichungssystem (15a) stellt die auf Mehrleitersysteme verallgemeinerte 
Form der iiblichen Leitungsgleichungen dar. Dabei ist der Punkt z = 0 auf 
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den Leitungsanfang bezogen. Legt man statt dessen z = 0 an das Leitungs- 
ende, wie es bei Leitungsproblemen oft vorteilhaft ist, so ist in (15a) das 
+-Zeichen an Stelle des —-Zeichens einzusetzen. 


2.3. Die aligemeinen Randbedingungen eines n-Leiter-S ystems 


In den Gleichungssystemen (15) und (15a) treten jeweils 2 Konstanten 
U’, U" bzw. U,, J, auf, welche aus den im speziellen Problem vorliegenden 
Randbedingungen abzuleiten sind. Physikalisch ist das Leitungsproblem durch 
die Abschlusswiderstande der einzelnen Leitungen sowie durch aussere Span- 
nungsquellen bestimmt, wie dies Figur 2 zeigt, welche ein Zweileitersystem mit 


She ‘ J}. 
a FZ 
Fig. 2 
Schema zur Aufstellung der Randbedingungen eines Zweileitersystems mit einem Null-Leiter. 


| Rie 
=t 


Null-Leiter, also das einfachste Dreileitersystem darstellt. Unter sinngemasser 
Erweiterung der Bezeichnungsweise der Figur 2 auf ein n-Leiter-System er- 
geben sich die folgenden Beziehungen zwischen den Konstanten der Gleichungs- 
systeme (15), (15a) und den vorgegebenen Werten der Abschlusswiderstande R 
und Spannungsquellen V 


U=-RAtVv, U=RiN+ V,, 


worin die Widerstandswerte R,;, Rj, zu Diagonalmatrizen und die Spannungs- 
quellen Vj; und V;, zu Kolonnenmatrizen zusammengefasst sind. 

Setzt man nun die aus (15) folgenden Werte ein, so erhalt man schliesslich 
die allgemeinen Randbedingungsgleichungen: 


(1aeR, G) Use Ry. GU" = Ve. 
ei”! (1 — R, G) U'+ e171 + R, G)U" = V,, 


(16) 


worin 1 die Einheitsmatrix bedeutet und zur Abkiirzung G = v ¢ gesetzt ist. 
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In ahnlicher Weise ergeben sich fiir das Gleichungssystem (15a) die Rand- 


bedingungsgleichungen: 
U, + Ry Jo= Vo, 


U, cosyl+7R,Gsinyl—j7 J, G-sin yl + R,cosyl = V,j. 


(16a) 


3. Anwendung auf das Dreileitersystem 


Verlaufen zwei parallele Leiter innerhalb eines Aussenmantels, so haben 
wir es mit einem allgemeinen Dreileitersystem zu tun, das nur im speziellen Fall 
geometrischer Symmetrie (gleiche Leiterquerschnitte, symmetrische Lage inner- 
halb des Aussenmantels) durch den Ansatz einer sogenannten Gegentaktwelle 
beschrieben werden kann. Ist eine solche Leitung elektrisch unsymmetrisch 
abgeschlossen, so lasst sich bekanntlich der Gegentaktwelle eine Gleichtakt- 
welle iiberlagern [3], und die Summe von Gleichtakt- und Gegentaktstr6men 
und -spannungen bildet ein System, das auch ohne Zuhilfenahme der oben 
abgeleiteten Beziehungen eine strenge Behandlung aller Probleme erlaubt. 
Diese Vereinfachung wird aber hinfallig, wenn die geometrische Symmetrie 
gestort ist, das heisst wenn die beiden Leiter verschiedene Querschnitte auf- 
weisen oder ihre Lage innerhalb des Aussenleiters unsymmetrisch ist. Es ist 
dann unter allen Umstinden zur Bestimmung der Strom- und Spannungsver- 
teilung auf den Leitern ein vollstandiges Gleichungssystem (15) oder (15a) auf- 
zustellen, dessen Lésung mit Hilfe der speziellen Randbedingungen zu er- 
folgen hat. 

Fiir den Fall eines Dreileitersystems nehmen die Gleichungen (15) die fol- 
gende Form an (mit G = v e) 


Jy) Gy, Gy» —Gyy 


—G, 2 | Corare 
BLAU] 
Jn |_| Gin Gre —Gin —Gee| | Wer 17 
U, 1 0) 1 ) | Uy" tive (17) 
: ae ys 1 0 eae U,! tiv | 


Auf ahnliche Weise ergibt sich fiir das System (15a) unter Anwendung der 
Regel zur Bildung inverser Matrizen: 


aI ki a COS z 0 —7Gysinyz —7G,siny z Tio 
Je 0 cosy z —7G,.sinyz —jGysinyz eb 

. ea . . G 9 . 
U, | =| -7 4 Sings +4 bs siny 2 cosy 2 0 | ee 

cr : _G ; 
U; +4 a sinyz —j See siny 2 0 cosy z bap 

; (Gy Gy» 
mit d= pox cee U (C41 Cog — Chg) 
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Je nach Problemstellung wird man eines dieser beiden Gleichungssysteme 
bevorzugen, und wir werden im folgenden fiir beide je ein Beispiel anfihren. 


4. Der Richtungskoppler 
4.1. Problemstellung und allgemeine Ausdriicke 


Unter Richtungskopplung verstehen wir das Zusammenwirken zweier Lei- 
tersysteme in der Weise, dass ein einseitig gerichteter Energiefluss im einen 
System eine verschiedenartige Erregung der beiden Enden des anderen Systems 
zur Folge hat. In der gewohnlichen Leitungstechnik zerlegt man die Kopplung 
zweier Leitungssysteme [4] in die magnetischen und elektrischen Kopplungs- 
anteile und stellt bekanntlich fest, dass sich diese beiden Anteile bei ebenen 
Problemen in der einen Richtung gegenseitig kompensieren. Als praktische 
Folge wird auf diese Weise zum Beispiel das Nebensprechen zweier paralleler 
Telephonfreileitungen in der Vorwartsrichtung sehr klein, wahrend das so- 
genannte Gegen-Nebensprechen erheblich gréssere Werte annimmt. Diese langst 
bekannten Tatsachen lassen sich in der Kurzwellentechnik zur Realisierung 
eines Richtungskopplers benutzen, der zum Beispiel eine direkte Ablesung des 
Reflexionsfaktors von Antennenzuleitungen erlaubt. Ohne auf den Mechanis- 
mus der Kopplung im einzelnen naher einzugehen, lassen sich simtliche Fragen 
eines solchen Richtungskopplers aus dem Gleichungssystem (15) unter Zuhilfe- 
nahme geeigneter Randbedingungen unmittelbar ableiten. 

Legen wir ein Dreileitersystem nach Figur 2 zugrunde, so haben wir folgende 
Fragen zu beantworten: 

a) Vorausgesetzt sei eine rein fortschreitende Welle der Spannungsamplitude A 
auf Leitung 1. Welches sind die Bedingungen, dass die Leitung 2 bei sym- 
metrischem Abschluss ihrer beiden Enden nur am einen Ende eine zu A 
proportionale Spannung aufweist, wahrend das andere Ende strom- und 
spannungsfrei bleibt ? (Bedingungen der Richtwirkung und Kopplung.) 

b) Unter welchen Bedingungen (Abschlussimpedanzen, Kopplung und Fre- 
quenz) kann iiberhaupt eine rein fortschreitende Welle auf Leitung 1 laufen ? 
(Anpassungsbedingung.) 

c) Wie lassen sich die Bedingungen a und b zur Verwirklichung eines Rich- 
tungskopplers vereinigen ? 

Da die Leitung 1 allein und nur am Anfang (z = 0) erregt werden soll, so ergibt 

sich fiir die Quellspannungen V, und V, nach den Bezeichnungen von Figur 2 


Vie ls) Vj=0: 


Setzen wir diese Werte in (16) ein, so erhalten wir die Bestimmungsglei- 
chungen fiir die Konstanten WOU; U,' in Matrizenform: 
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el — Ryu Gy) eM (L+ RG) -e 7” Ry Gye eV! Ry Gye 
=e"! Ris Gi ere Rs Gy, (1 — Ri» Go.) Opi cida 6) + Ryz Gy») 
Ros Gy. — Roe Gy» as: Ros Go. 1 — Roe Gyo 
1+ Ro Gy 1 — Ro Gy Rou Gis — Ro Gye 
(UL 0> 
U;' 0 
x] ae an (18) 
OTS) Ie 


Aus diesen Gleichungen lassen sich samtliche Eigenschaften der betrachteten 
Anordnung entnehmen. 


4.2. Richtwirkung 


Mit Richtwirkung bezeichnen wir die Méglichkeit, durch eine Messung an 
der Nebenleitung 2 die in der Hauptleitung 1 fortschreitende von der riick- 
kehrenden Welle zu unterscheiden. Als Mass der Richtwirkung dient das kom- 
plexe Verhaltnis der Spannungen an Anfang und Ende der Hilfsleitung 2: 


Ure 
=—H, 19 
as (19) 
Bei Einfithrung der Reflexionskoeffizienten auf den beiden Leitungen: 
in Uy 
Rea. Gee (20) 


lautet unsere Aufgabe, J} als Funktion von I, zu bestimmen. Im Dreileiter- 
system gilt nach (15) allgemein 
Use —= (Be ev! f Ue Use = U; + ie 


Fihren wir nun J} ein und kiirzen UZ, so wird: 


~? 


am aan rier 
<i paca 21) 


Nun ist J} in Funktion der eigentlich zu messenden Grésse I) zu berechnen. 
Dazu gehen wir von (18) aus, fithren Uy = I, Uj ein und setzen ferner die 
beiden Abschlusswiderstande Rog und Rj, einander gleich: Roe = Rip = Ry. 
Nach langerer Zwischenrechnung, die hier libergangen sei, erhalt man dann: 


pe Ae Ba Cua) (1 19! igre) 1 RG =) ent? 
2 = . - . 
pho Ba Gua) (Lo Th) oR ad RG vier ar get 


(22) 
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Zur Anpassung der zweiten Leitung schreiben wir vor: 


1 il 1 
Goo U Coe UV (Cop + Cyp) ’ 


Ry 


(23) 


das heisst, R, soll gleich dem Wellenwiderstand der Leitung 2 sein. Dann 
ergibt (22): 


—2jyl 
Di Cpe 


R= (24) 
Gleichung (24) in (21) eingefiihrt, ergibt schliesslich: 
pee hy Fale (25) 
Upe we pe a) 


Der Quotient |v| der beiden Spannungsbetrage U,, und Up, stellt also bei den 
gegebenen Randbedingungen unabhangig von der Lange der Leitung 2 ein di- 
rektes Mass des Reflexionskoeffizienten |J{| der Hauptleitung dar. 

Eine Welle, die in der Hauptleitung vom Anfang z = 0 gegen das Ende 
z= lauft, ruft dabei in der Nebenleitung Strom und Spannung am Anfang 
z = 0 hervor, wahrend das Ende z = / strom- und spannungsfrei bleibt. 


4.3. Berechnung des Kopplungsgrades 


Zu berechnen ist das Verhaltnis der Amplituden in der Haupt- und Neben- 
leitung: 
Uns 


p= : 
Ups 


(26) 


Es geniigt, dieses Verhaltnis fiir den Fall zu kennen, dass T, = 0 ist, und wir 
bezeichnen dies mit fy. 


Nun ist: 
Ups = U, + UZ = UZ (1 + Ji), 


Uy = Uj + U =U; (1+ Qh), 
und man erhilt fiir ky unter Beriicksichtigung der Anpassungsbedingung (23) 


Gi 
D Giro 


ko= (l= e- 77"). 
Ersetzen wir die Induktionskoeffizienten durch die Direktkapazitaten 


1 il 
as Gye = Cy, | Gon = Cao + Cys; 


ee 


282 Hans J. VON BAEYER und RONALD KNECHTLI ZAMP * 


so erhalten wir schliesslich fiir den Betrag von ky 


[to] = ct — |siny Z|. (27) 


Cre + Coo 


Man sieht daraus, dass eine optimale Kopplung bei //A = 1/4, 3/4, 5/4, ... zu 
erhalten ist und dass diese Kopplung in weiten Grenzen unempfindlich gegen 
Frequenzanderungen ist. 


4.4. Anpassung der Leitung 7 


Bis jetzt haben wir J} als gegeben vorausgesetzt. Wir untersuchen nun den 
Zusammenhang zwischen J} und den Randbedingungen; dazu greifen wir auf 
Gleichung (18) zuriick und erhalten nach einigen Zwischenrechnungen 


re io —2jyl 
in 1+FyiGy, Gee (28) 
1 Ry Gy. “a 12 (1h ‘g-2ir)} 
1 + Rr Gy 2 Goo ) 
mit 
o¢y1 11 Gi1—-1 
eg tet A ae 
of iat Ry, Gyyt+ 1° be 


Im Ausdruck (29) erkennen wir den Reflexionskoeffizienten einer gewohnlichen — 
Leitung (Zweileitersystem) ; und in (28) erscheint [4 einfach als ein in Funktion — 
der Kopplungskapazitat C,, = —(1/v) G,, korrigierter Wert von J],. Um zu — 
erreichen, dass die Einfiihrung einer Leitung, auf welcher J} zu messen ist, — 


dieses Jj nicht beeinflusst, muss der genannte Korrekturfaktor klein genug 
sein, das heisst, es muss im gefahrlichsten Fall der Anpassung der Haupt- 
leitung (R,, = 1/G,,) 


Unter dieser Voraussetzung ist die Anpassung der Leitung 1 und entsprechend 
die Messung von J} praktisch frequenzunabhangig. Allein die Kopplung | Ro 
ist frequenzabhangig, aber, wie wir gesehen haben, nicht in kritischem Masse. 


4.5. Verwirklichungsméglichkeiten 


Ein Dreileitersystem stellt nach den vorangehenden Uberlegungen dann 
einen weitgehend frequenzunabhangigen Richtungskoppler dar, wenn: 

Leitung 1 beliebig abgeschlossen und einseitig erregt ist, 

Leitung 2 beidseitig mit R, = 1/G,, abgeschlossen, das heisst angepasst ist, 


‘gh de toe 


12 
= ——* if 
GAG oe 


und / ungefahr = 4/4 oder ein ungeradzahliges Vielfaches davon gewahlt wird. 


~ 
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Zur Verwirklichung kann zum Beispiel innerhalb einer koaxialen Leitung 
eine zirka (A/4)-lange Hilfsleitung nahe der Aussenwand angeordnet, an ihren 
beiden Enden durch den Aussenleiter hindurchgefiihrt und reflexionsfrei ab- 
geschlossen werden. Diese Form wurde bereits kurz in der Literatur beschrie- 
ben [5], wobei jedoch die theoretische Begriindung wenig stichhaltig und zur 
Ableitung quantitativer Beziehungen ungeeignet war. Eine andere Realisie- 
rungsmoglichkeit stellt der sogenannte Schlitzkoppler dar, bei welchem zwei 
parallel gefiihrte koaxiale Leitungen durch einen schmalen Langsschlitz mit- 
einander gekoppelt werden. Dieses System wurde von KADEN [6] theoretisch 
behandelt, und es zeigte sich, dass, abgesehen von speziellen Problemen der 


| 
: | 
Zeal z=0 
Fig. 3 
Leitungsersatzbild der gefalteten Dipolantenne. 


Schlitzerregung, genau die gleichen Wirkungen erzielbar sind. Nicht zu ver- 
wechseln ist dieser Schlitzkoppler jedoch mit-den bei Hohlleitersystemen be- 
kannten Langschlitzkopplern [7], [8], bei welchen sowohl der Richtungssinn 
der Kopplung als auch der Kopplungsgrad vom Ort und der Breite der Schlitze 
sowie von der Wellenart im Innern der Hohlleiter abhangig sind. Dies rihrt 
daher, dass in Hohlleitern elektrische und magnetische Felder réumlich von- 
einander getrennt sind, so dass die Kopplung der Leitungen zumeist tiberwie- 
gend elektrisch oder iiberwiegend magnetisch ist, was zu wesentlich anderen 
Richtwirkungen fiihrt als die vorliegend beschriebene elektromagnetische Lei- 
tungskopplung. 


5. Die gefaltete Dipolantenne 


Ein weiteres Beispiel fiir die Anwendbarkeit der Mehrleitertheorie in der 
Ultrakurzwellentechnik stellt die gefaltete Dipolantenne (folded dipole) mit be- 
liebigem Durchmesserverhaltnis der beiden Leiter dar, deren Stromverteilung 
sich in erster Naherung leicht aus den Gleichungen (15a) ergibt. Als erste 
Naherung verstehen wir dabei, wie stets in der Antennentheorie, die approxi- 
mative Bestimmung einer Stromverteilung durch Anwendung der Leitungs- 
gleichungen unter Vernachlassigung der Abstrahlungsverluste. Die Erfahrung 
zeigt dabei, dass diese Losung ein recht genaues Bild der grundsatzlichen Strah- 
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lungseigenschaften ergibt. Zur Erfassung des Impedanzverlaufs muss selbst- 
verstandlich die Abstrahlung beriicksichtigt werden, jedoch auch hierfiir lassen 
sich approximative Angaben machen, die zu befriedigenden Ergebnissen fiihren. 

Als einfachstes Ersatzbild einer gefalteten Dipolantenne geniigt es, nur die 
eine Dipolhalfte nachzubilden, und man erhalt die Leiteranordnung Figur 3. 

Zur Berechnung ist die Verwendung des Gleichungssystems (15a) vorteil- 
haft, das fiir drei Leiter (Leiter 1, Leiter 2 und Erde 0) die Form (17a) an- 
nimmt. Die Randbedingungen sind der Figur 3 zu entnehmen. 


Upi = Ue, Jor = —Soz: Uj2=0, Us=Vv. 


Statt der letzten Bedingung fiihrt man vorteilhaft die Grésse J = J, an der 
Stelle z = / als Bezugswert ein, da dieser Strom auch bei beliebigen Langen / 
niemals Null werden kann, wahrend die ungedémpfte Resonanz des Systems 
bei bestimmten Langen zu einem Kurzschluss am Leitungsende mit U,, = 0 
fihren kann. Unter Beriicksichtigung dieser Randbedingungen und unter Ver- 
wendung der Direktkapazitaten nimmt das Gleichungssystem (17a) folgende 
Form an: ; 
a 


1+ Acos?y/ 


\{— (A + 1) cosyl cosy z + a siny J siny z\, 
20 


1+ Acos*yl 
a (30) 

Sot leak. 
1+Acos?yl 


- | == sin y L.cosy ob eae Cus 1 sin 2| 
mire Cig PERN ae 


20 


Ree, 


| | 
h={ J | {(4 + 1) cosy Dcosy «+ siny Usiny Z}, 
| 
| 
| 
| 


a Ae hie aa 
ix a | {siny / cosy z — cosy /siny z} 


A=Ci (= <-): 


Hiermit sind sémtliche Eigenschaften des in Figur 3 zugrunde gelegten Leiter- | 
systems bestimmt. Fiir die Strahlung kommt von der Stromverteilung auf den 
beiden Leitern 1 und 2 nur derjenige Teil in Frage, der ein beide Leiter um- 
schliessendes Magnetfeld besitzt. Dies ist einfach die Summe der beiden Stréme 
J; und J,: " 

At jp af (1 + at) (Peas) siny z. (31) 
Hieraus sieht man, dass sich die gefaltete Antenne bis auf einen von der An- 
tennenlange / abhangigen Faktor genau gleich wie eine normale Stabantenne 
verhalt, namlich sinusférmige Stromverteilung aufweist. Bei yl =2/2 (Ab- 
stimmung auf eine Viertelwellenlange) ist dieser Faktor gleich 1, und der ge- 
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meinsame Strom (J, + J) ist um den Faktor (1+ Cy9/Cy») grdsser als der Speise- 
strom des Leiters 1 am Orte z =/. Bei gleicher abgestrahlter Leistung muss 
somit der durch Strahlung verursachte Anteil der Antennenimpedanz am Spei- 
sungsort bei der gefalteten Dipolantenne um den bekannten Betrag (1+ Co9/Cy9)? 
grdsser sein als bei der Stabantenne [9]. Auch der Blindwiderstandsgang der 
Eingangsimpedanz ergibt sich aus (30) durch die Bildung des Verhiiltnisses: 


Cc 
if 20 
Sani ag, (32) 
Ji UCyo tgyl—v Cy (A+ 1)ctgyl ° 
Fir A = 0, das heisst Cj, = 0 und Cy) = Coo, geht dies iiber in 

U; 1 

Se ee te ala 

eee tuo (33) 


Dieses Ergebnis darf aber fiir den Fall der strahlenden Antenne nicht unmittelbar 
ausgewertet werden, da bei seiner Ableitung der ohmsche Anteil vernachlassigt 
wurde. Eine approximative Bestimmung des Impedanzverlaufs wird aber durch 
Aufstellung eines Ersatzbildes méglich, welches in der Parallelschaltung einer 
ungedampften, kurzgeschlossenen und einer gedampften, offenen (A/4)-Leitung 
besteht. Dies wurde in der Literatur [9] bereits beschrieben, so dass sich ein 
weiteres Eingehen darauf an dieser Stelle eriibrigt. Es sei nur darauf hinge- 
wiesen, dass Gleichung (32) fiir den Fall verschwindender Kopplung zwischen 
den beiden Leitern (A = 0) bereits deutlich den Impedanzverlauf zeigt, der . 
zur Aufstellung des erwahnten Ersatzbildes fiihrt. 

Die vorliegende Untersuchung entstand im Mikrowellenlaboratorium der 
Aktiengesellschaft Brown, Boveri & Cie. in Baden als Resultat verschiedener, 
iiber mehrere Jahre verteilter Studien. Eine experimentelle Auswertung der 
mitgeteilten Methoden fiihrte dabei zu einer Anzahl im praktischen Gebrauch 
bewahrter Konstruktionen. 
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Summary 


The paper deals with problems of wave-propagation on multiconductor 
transmission lines at very high frequencies. The transmission line equations for 
current and voltage on each conductor are given in a generalized form using 
matrix symbols. Two examples show how these equations are to be applied to 


practical problems and how proper boundary conditions are to be introduced. — 
The first example concerns the theory of a transmission line directional coupler, — 


the second example a simplified theory of the folded dipole antenna. 


(Eingegangen: 18. Februar 1952.) 


Note on the Fundamental Solution of w y,,+ y,,,=0 


By JuLian D. Core!), Pasadena (California) 


1. Introduction 


The equation 
W Vov a Yow = 0 (1) 


occurs in the gasdynamical theory of transonic flow. It is of mixed type, 
elliptic for w > 0, hyperbolic for w < 0. In this note a brief study is made of 
the fundamental solution of (1) using ideas related to the physical applications. 
Hence no attempt will be made at a precise discussion of conditions on func- 
tions, etc. except where necessary. The basic problem here is to obtain a re- 
presentation of the fundamental singular solution of (1) with a singularity in 
the elliptic half plane w > 0. The nature of the singularity and the continuation 
of the solution into the hyperbolic half plane w < 0 are considered. It is in 
connection with this continuation that the physical problems provide some 
insight. 

There is already in existence a considerable mathematical literature on (1) 
starting with the historic memoir of Tricom1[1]2). Most of the references on 


1) Guggenheim Aeronautical Laboratory, California Institute of Technology. 
*) Numbers in brackets refer to References, page 297. 
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the subject are included in [1] which contains a bibliography, and several others 
are appended here ([7] to [10]). In addition various authors have developed 
results for special physical problems. However, it is believed that the methods 
used here, although by now standard in the theory of partial differential equa- 
tions, are sufficiently different from those appearing in the literature to be of 
interest. Some new results are also presented. 

Some of the results appearing in this note have already been used in a 
previous paper [2] and some are used in a paper to appear soon [3]. In both 
applications doublets formed from the fundamental singularity were of interest. 
These doublet singularities represent, in the hodograph, the flow at infinity. 
It will appear in the analysis that the solutions formed from the fundamental 
singularity should be useful when the hodograph is one-sheeted at infinity. 
The case of branched singularities is reserved for [4] and future work. 

The author thanks T. Y. Wu for his generous assistance on many points. 


2. Definition and Representation of the Fundamental Solution 


Consider 
W Voy a Ae Q(w, v) (2) 
in the domain 
Etoo UW =< co, — oo < P< Go. 


A fundamental solution of (1), T(w, v — v,; w,) is defined as a singular solution 
of (1) such that 


y(w, 0) = ff Tw, v — 043 0) Ow 2) di, doy 3) 


is a solution of (2). It follows from the invariance of (1) under translation with 
respect to v that T is a function only of the difference v — v,. A definition 
equivalent to (3) is that 7(w, v; w,) is a solution of 


W Von + Vow = —d(v) O(w ie 4) S (4) 
It is assumed that w, > 0 as this is of most interest for applications. w, = 0 
will be treated as a limiting case. 
An integral representation of the fundamental solution T(w, v; w,) can be 
obtained by Fourier analysis of (4) with respect to v. Let 


co 


ee fei ye V5; w,) dv (5) 
7 


y(w, B; w,) a 
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so that 


y(we, v5 10) = | oF Hee, Bi te) vB (6 


Applying the transform (5) to (4) yields 


— oo 


1 


Le cei toma 
Vow — Pw y Vin 


6(w — w) . (7) 


The Fourier transform of the fundamental solution T is thus the fundamental 
solution of the ordinary differential equation (7) (divided by 2 z). This solution 
for ordinary differential equations is discussed in [5] (vol. I, p. 302). A solution 
of (7) is thus defined by the conditions 


(a) y is a continuous function of w; 


(b) y is a solution of y,,,, — 8? wy = 0, except 


at w= w,, and y vanishes at +00; t (8) 
(c) jump condition 

; dy _ (@ oe ©: 

‘at (as jeer (a ee - V2 : 


The solutions of (7) for w > 0 are linear combinations of the modified Bessel 
functions 


Bs Ky) (\B| Sie At Ty 1/3 (|B 2) (9) 


where z = 2 w*/?/3. The notation is that of [6]. The conditions (8) do not in 
this case determine the solution uniquely. However, the arbitrariness which 


occurs may be exhibited and the Fourier transform of the fundamental solution 
can be written out as 


Cy (=) Ky)y (|B| 2) (z > %) 


| Cy ( =)" K,, (|B| 2) + 3 (5)" (10) 
x [Lis (|B | 2) I 1/3 (|B| a Ty) (|B | 2) Tig (\B| 2)]. 
(zie) 


V2x T(w, B;.w;) = 


Here C, = C,(z,, 8) is arbitrary to the extent that the inverse transforms exist. 
The term in the solution multiplied by C, is regular while the bracket in (10) 
represents the singular part. Thus although many fundamental solutions exist 
the singularity is the same in all of them. C,(z,, 8) can be chosen so that the 
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fundamental solution satisfies some supplementary conditions. These conditions 
are concerned with the behavior of the solution in the hyperbolic region w < 0. 


3. Continuation into w <0, Supplementary Condition 


The continuation of (10) into the hyperbolic region w < 0 is provided by 


eT (pl) =o Jay (PLO: | 
PPT. pl Bz) = 0® Foie (BLO: 


2° Kyy (\B\ 2) =o" a [F 13 (BLS) + Says (1B) 9] | 


(11) 


where 


D i ‘ Se i/2 
paw, Fa PF, w=—u, or 2S te**. 


(10) becomes 


oe TC Bs2,) = ear ier. c, — 2° 1, (\B| a)| Tan (BIS 


(BY (a)? G88 Law 81d] Jin 1819 - 
02 


It is convenient to split C, into its odd and even parts with respect to B. Since T 
is real let 

Cyler B) =[Cale, B) it Coley BZ) (13) 
where 


Co(%1, B) = Colt, — B), Cs(t, B) = —C3(%1, —8) - 


Using the inversion formula (6), T is represented by the following integrals 
for w < 0 


T(C,v;%) = = & oe iB [MG.- 2" Ly) (B 21) Ta 1)s (BC). 


It 3 


+ (CG. ai Typ (B 21) Sips (B ¢)} cosB v dp (14) 


3 


4+ = (- a jx | Clam (B 6) + Sis (8 £)] sin B v ap. 


The type of supplementary condition to be introduced is familiar from the 
theory of purely hyperbolic equations. For hyperbolic problems the direction 
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of propagation of disturbances along the characteristics, or the time-like direc- 
tion, must be prescribed in addition to the differential equations. Using this 
idea, regions of influence and domains of dependence can be specified. For the 
equation of mixed type (1) a similar requirement on the behavior of the solution 
in the hyperbolic region makes the fundamental solution T unique. 

A physically significant condition is suggested as follows. The lines y 
= constant may be thought of as corresponding to the streamlines of the flow. 
The two families of characteristics of (1) are 


Tee = us? — y + € = constant, 
(15) 
VP. . 


2 fc P 
os u*? — y — € = constant, 


and these are the hodograph of the Mach waves. The time-like direction is as- 
signed as the direction of flow along the streamline. Thus the direction of prop- 
agation of disturbance along the char- 
acteristics or Mach waves emanating 
from a point on the streamline is 
assigned. One important type of flow 
is locally vortex-like. In this flow one 
of the Mach waves leaving the stream- 
line belongs to a family of expansion 
waves and the other to a family of 
compression waves. (The other types 
of flow which can occur locally are 
source-like, both characteristics are 
expansion or both compression, or 
characteristic-like where the stream- 
line fallsona characteristic.) A typical 
Fig. 1 hodograph picture of vortex-like flow 
is shown in figure 1. Hence, a direc- 
tion of disturbances is assigned to the characteristics. On I_ disturbances run 
from u = co toward uv =0 and on I’, from « = 0 toward u = oo. This state- 
ment of direction of propagation along the characteristics is the necessary sup- 
plementary condition. 

The condition may be formulated on the behavior of y as u > oo. If the 
values of y are caused by a singularity in the elliptic half-plane w >0 no 
disturbance should exist on the characteristics "_ as u > oo. That is, the values 
of y should depend only on the values of I’, as u > oo. It turns out that this 
condition can only be realized approximately (there is always a certain amount 
of reflection on J") and that it may be required 


y= lM tn (vt t) as $00. (16) 
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It may be observed that this condition introduces unsymmetry of the solution 
with respect to v. 
4. Representation of the Fundamental Solution 


The condition (16) suffices to determine the functions C,, C3 of (14). An 
asymptotic approximation to (14) as ¢ + oo is obtained by replacing the 
Bessel functions with the first terms in their asymptotic expansions 


pee |e cos(BE— 35), | 


Tig BO Va ep 003/80 — 55). | 


Then the integrand in (14) contains products of trigonometric functions which 
may be combined to contain terms like 


easoige-c) cusp eC )o “sing (7 —%)., Sing 10 4 C) 


If, according to (16), the coefficients of the functions of v — € are set equal to 
zero, C, and C; are determined as follows 


Cr @6) = [hp Ba +L, Ga), 86> (8) 


Cs (@ B) = — + > KiplB 4) - (B > 0) (19) 


Thus the following Fourier integral representations of the fundamental solu- 
tions are obtained by substituting the values of C, and Cz in (14), (13), and (10) 


| 


(17) 


T(e, 03) = pe (752)” [Kip B29 [Lin 8) + Lin BaD] 


xicos BU — = Kyyg (B 2) sin B of ap , (z > 2) 


=k (228) [Ky Bad (an 9) + Ln GA) 


x cospv — — Kyyy (B 2) sin B 0} ap , (gai) 


Die, v; 4) = = ( : 2 ay f Kap (B 2) {es (Bp (ee Jiy3 (8 ¢)] 


1 


75 Fs B 0) + Sin B O1sinb ol dB. (> 0) | 
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~ eal 


It should be noted that if doublets or higher singularities are formed from | 


T of (20) by differentiation with respect to v or z,, these doublets also satisfy 


the equation (1) (except at the singularity) and have the same behavior as T 


in regard to propagation of disturbances along the characteristics. 


5. Singularity at w, = 0 


As a check on the ideas of the previous paragraph, the limiting case w, = 0 
can be considered. According to the theory of characteristics, if there is a 
discontinuity (in derivatives, say) across a characteristic at one point then a 
discontinuity persists all along the characteristic. The same can be said for a 
singularity of the solution. Thus if there is a singularity at « = 0, v = 0 there 
may be, in general, singularities on both the characteristics v -+ € = 0, v—£=0. 
However, since it has been required that disturbances propagate away from 
the sonic line u = 0 only on the characteristic J",: v + € = 0 it should be ex- 
pected that the fundamental solution (20) has no singularity on [_: v — £ = 0, 
but does have one on J",: v + ¢ = 0. This result is verified by evaluating the 
integrals in (20) which represent T. 

If z, > 0 in the integrals in (20) valid for £ > 0, (20) becomes 


ri o.= 69 = A (BP 
a 
: | {Fa (8 5) — Jijs (B £)] cosB v (21) 


— Fp Ua 82) + Ji (BO) sino} 2A 


The integrals occurring in (21) can all be evaluated in terms of hypergeometric 


functions and the formula is given on page 401 in [6]. Using some transforma- | 


tions of hypergeometric functions the results are tabulated below in a form 
exhibiting the singularity 


& | | 0 (0<¢< |v}) | 
fc0sB Jy (BO) BM ap = : } (22) 
0 | A C18 (4 aa ayes (0 <. |v| a ¢) 
Pe | A V3 fils ar £2) 176 (OC ||) 
| cosB Ty (BS) B- ap = (23) 


A2 Tae pais O< |v] <<) 


> 


' 


-Vol. III, 1952 Note on the Fundamental Solution of w y,,,, + Vy» = 0 293 


si 1/8 
A2 C118 amare a ° 
a © caw ag |x PUb ee Be a aa) OSE < Led 
Sf sinBo TB BP MB=) (24) 
) 
P AYV3 E18 (f2 — y2) 1/6 B ssw 
Dh a v? 
KF(L >, G 1-H) O< [1 <d) 
fos) A sign v Meera 
Bsiopo Ioan BIA aB=_ POO” pes | (25) 
0 0 (Or yu |<.) | 
where 
1 
A r(=) 2218 


ce Vx 2418” Bs r(S) 


The integrals (22) to (25) all have 1/6 power singularities on v + ¢ = 0. How- 
ever, the linear combination of these solutions which forms the fundamental 
solution T has no singularity on v = ¢. Substituting (22) to (25) into (21) yields 


31/6 22/8 ees r2 
Bix iagilh 5 1a) | 
i 
re ) 
1 1 6 Isr ea’ 
+ Sie ° p58 yi * z i aie an (0< ¢< |») 
(5) 
Z(G, 0) = { af ee i (26) 
sel v oe 
ang peF(L Ze 1a) 
1 
r=) | 
il il 6 1 — sign v) 
33 35/6 © 2518 x 1l2 * nN ro y2)ile (O< |2| < 6) 
"(3) 


The discontinuity in (26) at v = 0 is only apparent as can easily be seen by 
using the transition formula for the hypergeometric function 


oe 
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6. Nature of the Singularity for w, > 0 


It is also of interest to exhibit the singularity of T(z, v; z,) given by (20) 
at z= 2,, v = 0. It is known from the general theory of elliptic equations that 
the fundamental solution has a logarithmic singularity at z = z,, v = 0 and 
the exact nature of that singularity can be worked out in this case. 


Lee 
1 2 1/3 = 
2. (458)" [inno 
0 
xX {Lijjg (B 21) + Lig (8 21)} cosB v dB pheucca ma) 
Tee { 1/3 1 BS 1 } (28) 
iL 2 1/3 
te (AEE) fe 
0 
x {ia (Bz) + Tig (B z)} cos B v dB (z< 4) 
1 2 1/: an ™ ; 
Tle, 034) = — ae (—54)" |Z Ky BA) Kin (Bx) sinB v dB (29) 
0 
so that 
T = 7, J; (30) 


It can be seen from (29) that 7, = 0 if v= 0 so that, for z, > 0, all of the 
singularity of T is contained in 7,. The purpose of adding T, to T, is to satisfy 
conditions as € - oo. In addition, as seen in the previous paragraph, 7, 
contributes a certain type of singularity for the limiting case z, = 0. 

The integrals representing 7, in (28) become equal to known integrals if 
the integration in (28) is regarded as taking place in the complex plane and the 
path of integration is deformed to the imaginary axis. If this is done, (28) 
becomes a Laplace integral 


T,(z, 0; 2%) => REN 


Le, °) 


3 | (31) 
ne yh ell ST sin (2) Tapp (e%) + iyo (#2) Juyy (0 41) dn | 


0 


a representation valid for all z. The integrals in (31) are special cases of a formula 
in [6], p. 389, namely 


f Re (a+ib+ic)>0 
1 2 2 2 
[ ene J, (6t) f, (ct) dé = —__ Q,_,, (75), ; 
é a Vb c, Rey>—+ 32) 


— 


Vol. III, 1952 Note on the Fundamental Solution of w y,,. + Vy» = 9 295 


where Q,, is the Legendre function of the second kind. Q, is expressible as a 
hypergeometric function 


r Vx I'(n + 1) 1 mt+ti =n 3 il 
Qn(Z) =~ 9n+1 * T(n + 3/2) * 7n+1 F( 2 2D +1jn+ 555 i (33) 
Applying (32) and (33) to (31) we have 


T,(2, 0, 4) = T34(2, ¥, 4) + T-(@, 2, %) 


with 
1 / r 
Dee NES a. 
TENS, 0,85) = (4) Je Ie Jat peul Ose ind erate 24 )naee 
0 
eee! (e\s 1 v2+ Z2-+ 28 
oe o4 a2,\;3 Gaye O56 a) (34) 
il 
ap ae ee 1 re) oy Te al 
4yn \3 (222,Z)Ue © B fe: Sal 
‘ iP ey 
where now 
we v2 + 2274 23 
- Oe, 
and 
| co 
I ff Ae NYE ~ 
Doa(2, 0, 24) = (52) fet" Sys (2) ipa ) a 
0 
1 2\1/3 1 v2 + 224 ZF 
— Ag | CEA O10 ( Die ) (35) 


ee ts) 


uy _ (3)" Bas 
4 /x 3 CRA CPA de rs) 


Bic; ea oaets eal 
Ele LA atk 


From (34) and (35) it is evident that the nature of the singularity as v > 0, 
z > 2, is given by the singularity of the hypergeometric functions as Z > 1. 
To exhibit the singularity a formula is needed relating the hypergeometric 
function of 1/Z? to a function of 1 — 1/Z°. The necessary formula is 


F(a, bj a +b;0) = 3" om (1— 2)" [B, + log (1 — 0)] (36) 


where 
T(a+ 6) P(a+n) r(o+%) 


LL eel CHARA TR Pe 


Bn = ylatn) + pb+n)—2y +7). 
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(36) is an exceptional case for hypergeometric functions and can be derived 
by use of BARNES’ method ([11], p. 286)1). If (36) is used in (34) and (35) the 
complete expansion results 


1 1 1 


T,_(2, ¥, 2%) = — 40 (V2) 318 "“(u? + 22 4 ze” r(+z) r=) 
12 ale, 
1 r U \ 
oo r(n+ 5) (n+ 35 vt 2u2 (z24 22) + aaah (37) 
* ot Fin Fi) vit 208 (2 + af) + (2? + a8)? 
( Seta ea 
x {Bn + log v4 + 2 v2 (2? + 22) + (2? + 23)? 
where 


Br = va +) + olay +") -2y (+m), 


21/6 (z 2,)?8 1 
hae 22 v, 2) = Tipe eREy © (v2 + 22 4+ 22)5/6 i r 5 r a ii | 
(te) PGs) 
5 =) 
pitta) tig v4+ 2u3 (22 + 2%) + (22 — 29)? 
= T?(n + 1) v4 + 2u? (22+ 22) + (22+ 28)2 


| 
| (38) 
) \ | 


4 2y2 24 2%) 4 2__ »2\2 
x {Buy + log s+ u2"(2 22) + (2 ar} 


vi + 2? (22 + 2?) + (224 22 
where 


Bur = v(q> +m) + vay + 2) -2y (04m). 


The dominant term in the expansion near v = 0, z = 2, occurs for 7 = 0. Thus 


14 222, 2/3 
fs 1 + (pay) 
Oe tay x (V2) 38 (uo? fat agin 
ug v4 + 202 (22 + 2%) + (22 — 22)2 ne 1 (39) 


v4+ 202 (22 + 23) + (2? + 22)2 7 (V2) 31/8 
ie. Shoe (Bo) (2 2 2,)?/8 * 


“Vora ape + wag tp pare 


(39) shows that T has a logarithmic singularity of a rather complicated argu- 
ment (cf. [1], p. 68). (39) also shows how the logarithmic singularity disappears 


1) The formula was also supplied to me by Professor A. ERDELYt. It occurs in the (unpublished) 


Bateman Handbook of Special Functions as a special case of the formula for F(a,b;a+b+1;0) 
Where} == 0, 1012.05. 


pee = ame 


te 


Oe te 


ee 
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as 2, > 0 and a 1/6 power singularity remains. It should be remarked that if 
doublets or higher singularities are formed by differentiation of T with respect 
to v or w, the resulting expansion contains a logarithmic singularity in addition 
to higher singularities. 
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Zusammenfassung 


Die Grundlésung der Differentialgleichung w Vyv + Yow = 9 wird untersucht. 
Integraldarstellungen von Lésungen, die eine Sin gularitat im Bereich w > 0 besitzen 
und in einer besonderen Art in den Bereich w < 0 fortsetzbar sind, werden gege- 
ben. Die Singularitat fiir w = 0 kommt dabei als ein Spezialfall heraus. Die Rei- 
henentwicklung des singularen Teiles der Grundlosung wird aus der Integraldar- 
stellung abgeleitet, und der Charakter der logarithmischen Singularitat fir w > 0 
wird dadurch gezeigt. Diese Resultate finden ihre Anwendung in der Theorie der 
schallnahen Potentialstromung. 


(Received: April 7, 1952.) 
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Die automatische Kursregelung von Flugzeugen 


Von Rupo.F GLaus, Minneapolis, Minn.*) 


1. Einleitung 


Bei der Lésung von Regelproblemen wird im allgemeinen vorausgesetzt, 
dass das zu regelnde System durch lineare Differentialgleichungen mit kon- 
stanten Koeffizienten beschrieben werden kann. Diese Voraussetzung ist mei- 
stens geniigend genau erfiillt, solange man sich auf kleine Abweichungen von 
der Beharrungslage beschrankt. In diesem Aufsatz werden nur lineare Systeme 
behandelt. 

Beim Entwurf einer Regelung sind immer bestimmte Reglerelemente ge- 
geben. Die Lésung besteht dann darin, zusatzliche Elemente so zu bestimmen, 
dass zusammen mit den gegebenen eine vorgeschriebene Regelaufgabe erfiillt 
wird. Hierbei muss vorerst das dynamische und statische Verhalten der gege- 
benen Elemente bestimmt werden. Weiter ist zu entscheiden, welches die zu 
regelnden Grdssen sein sollen. 

Die Lésung von Regelproblemen wird iibersichtlich durch Aufzeichnen 
eines Blockdiagramms. Hierbei entsprechen die einzelnen Blécke den verschie- 
denen Regelelementen und enthalten ihre Eigenschaften. Falls die Anforde- 
rungen an die Regelung bekannt sind, kann mit der Berechnung der zusatz- 
lichen Systemkomponenten begonnen werden. Damit wird vorerst ein voll- 
standiges Blockdiagramm der Regelung erhalten, welches das statische und 
dynamische Verhalten des ganzen Systems beschreibt. Die letzte und wich- 
tigste Phase ist dann die Transformierung des Blockdiagramms in physikalische 
Komponenten. Die analytisch oder graphisch gefundenen Einstelldaten sollten 
bei sorgfaltig durchgefithrter Synthese nicht mehr als 30% von den im Ver- 
such gefundenen abweichen. Die Bestimmung der Einstelldaten von Steue- 
rungen auf rechnerischem Wege ist bei komplizierten Regelungen speziell vor- 
teilhaft. Nehmen wir zum Beispiel die vollautomatische Steuerung von Flug- 
zeugen, wo mehr als ein halbes Dutzend Parameter eingestellt werden miissen, 
die zum gréssten Teil nicht unabhangig voneinander sind, so ist einleuchtend, 
dass nur durch Probieren eine Optimaleinstellung sehr schwierig und zeitrau- 
bend, wenn nicht unméglich wird. * 


1) Minneapolis-Honeywell Regulator Company, Aeronautical Division. 
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2. Grundbegriffe 


Regelungen von der Art, wie sie hier untersucht werden sollen, kénnen 
durch ein Blockdiagramm dargestellt werden, wie aus Figur 1 ersichtlich ist. 

Dieses Blockdiagramm kann am besten an Hand eines Beispiels beschrieben 
werden. Es moge #,, der Eingang oder Sollwert zu unserem Regelsystem, die 
zeitliche Anderung eines Winkels darstellen. Diese zeitliche Winkelvariation #, 
sei der Ausgang eines Kommandogerates und bedeute das gewiinschte Azimut 
eines Geschiitzes. Die kontrollierte Regelgrésse 3, wird als Ausgang oder Ist- 
wert bezeichnet. Sie ist das tatsachliche Azimut des Geschiitzrohres. Die Aus- 
gangsgrosse #, ist durch die sogenannte Riickfiihrung zum Eingang zuriick- 
gefiihrt und wird dort mit dem Eingang verglichen. Damit erhalten wir ein 


REGLER | 2g |VERSTARKE SERVOS eoeeon wh 
MOTOR LAST 
Ky vig) Ks Gsip Kes Gea 


RUCK FUHRUNG 


Fig. 1 
Blockdiagramm eines Regelkreises. 


geschlossenes Regelsystem, auch Regelkreis genannt. ¢ sei die Differenz zwi- 
schen Eingang und Ausgang und werde kurz als Abweichung bezeichnet. Die 
Abweichung wirkt auf einen Regler und dessen Ausgang auf einen Verstarker. 
Da im allgemeinen die Leistungen der Eingangsgrésse #,, der zuriickgefiihrten 
Ausgangsgrésse %, sowie der Abweichung ¢ und des Reglerausganges @, sehr 
klein sind, brauchen wir einen Verstarker. Dieser wirkt auf einen Servomotor, 
der die Leistung zur Bewegung des Geschiitzrohres aufbringt. Das Geschttz 
wurde in diesem Beispiel die Regelstrecke genannt, wahrend als Last Storun- 
gen, die von aussen auf den Ausgang wirken, bezeichnet werden. Die Losung 
des Regelproblemes kénnte so umschrieben werden, dass das Azimut des 
Geschiitzrohres #,, dem kommandierten Azimut #, tiber einen vorgeschriebenen 
Frequenzbereich mit einer bestimmten Genauigkeit folgen muss, ungeachtet 
dusserer Stérungen. Somit stellen wir an die Giite unseres Regelsystems be- 
deutend héhere Anforderungen als zum Beispiel bei Drehzahl- oder Temperatur- 
regelproblemen, wo meistens der Sollwert ein Festwert ist und das System 
nur Lastainderungen unterworfen wird. 
Es ist jetzt fiir Figur 1: 
6H = Fey — Yaw ° i) 


Um nun den Ausgang als Funktion des Einganges fiir beliebige zeitliche Ande- 
rungen des Einganges beschreiben zu k6énnen, benutzen wir die Laplace- 
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Transformation [1], [2], [3]4). Damit kénnen wir schreiben: 


L [Pi] = a, (p) * (2) , 


Hierbei bedeutet L[Y,,)] die Laplace-Transformation von #,%. £ ist eine 


komplexe Variable: 


p=a+tw, (w= Kreisfrequenz in s-!, 1= aeey 
L[8,4] = 3, wy? L [é)] = Lp) (3) 
€) = Bep) — Faipy - (4) 


Durch die Laplace-Transformation kann zum Beispiel das Verhaltnis von Reg- 
lerausgang zu Reglereingang folgendermassen ausgedriickt werden: 


Kp Grip) = ae (P) - (5) 


é€ 


Hierbei wird Kp Gri») der Frequenzgang des Reglers genannt. Es ist dann 


oO 
ee (p) =* Kp Grip) ; K, Gy ip) a Ks Gs (p) * Krs Crs p) mart G »)- (6) 
K Gi) wird als der resultierende Frequenzgang des aufgeschnittenen Regel- 
kreises bezeichnet. Das Ausschreiben des Frequenzganges in der Form K Gp) 
hat seine Begriindung darin, dass der Frequenzgang in einen frequenzunab- 
hangigen, konstanten Teil K und einen frequenzabhangigen Gi) aufgespalten 
werden kann. K wird der Verstirkungsfaktor genannt. Mit Gleichung (1) bis 
(6) erhalten wir: 

u SoG; 

pr Mesiren dam: f (7) 


Gleichung (7) erlaubt, aus der Gangkurve K Gi») des aufgeschnittenen Regel- , 


kreises, den Ausgang als Funktion des Einganges fiir den geschlossenen Regel- 
kreis zu berechnen, bei Abwesenheit einer Last. Soll ®, ») fiir eine beliebige Ein- 
gangsfunktion #,, bestimmt werden, so kann folgendermassen vorgegangen 
werden: 

1. Bestimme die Laplace-Transformation von O. 

2. Berechne den Frequenzgang K Gy) » 

3. Mit Gleichung (7) bestimme (&,/0,) (p) , 


4. Schliesslich wird mittels der inversen Laplace-Transformation, die wir 
symbolisch als L~? schreiben, 


t)» 


1) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis auf Seite 308. 


— 
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KG 
Bigg hoe (p) 
2G Ee aie (8) 
erhalten. Falls der Eingang sinusférmig variiert wird, so ist der Ausgang eben- 
falls durch Gleichung (7) bestimmt. Nur ist 7 m an Stelle von # einzusetzen: 


KE Gj.) 


Baie) =e 1+ KG, 


o 


e(iw) (9) 


1m) 


(0,/0.) (tw) wird wieder als der resultierende Frequenzgang des aufgeschnit- 
tenen Regelkreises bezeichnet. 
Wenn der Eingang ein Stoss ist, erhalten wir: 


(10) 
und mit Gleichung (7) : 


cA @ 1 


Die rechte Seite von Gleichung (11) ist eine rationale Funktion und kann des- 
halb auch folgendermassen geschrieben werden: 


Ane eA zp i An ; 
a i er ae oe Oe 

Hierbei sind #,, fs, ..-, Py die negativen Wurzeln der Gleichung 
1+ KGS 0: (13) 


Die inverse Laplace-Transformation auf Gleichung (12) angewendet, ergibt: 
8, = Aot Ax EP A, et es Ae tn (14) 


Aus Gleichung (14) ist ersichtlich, dass der Ausgang mit zunehmender Zeit 
die Konstante A, annahert, falls alle Wurzeln fj, ..., fy entweder positiv reelle 
oder komplexe Zahlen mit positivem Realteil sind. 


3. Eigenschaften des Frequenzganges K Gio) 


Der resultierende Frequenzgang K G;;,,. wird erhalten durch Ersetzen von 
p durch 7 w in Gleichung (6): 


= 28 (iw). (15) 


Hierbei ist (#,/e) (i w) eine Vektorgrésse und kann beschrieben werden durch 
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das Amplitudenverhaltnis A von Ausgang tiber Abweichung: 


A =|"* Go) (16) 


und die Phasenverschiebung 


p= x (ia). (17) 


Das Zeichen < (#,/e) (¢ w) bedeutet den Winkel zwischen den Vektoren %,(;..) 
und é;;.,) . Die Gangkurve K G;;,,) ist in Figur 2 in Form eines Polardiagramms 
dargestellt, wobei die Frequenz w als Kotierung der Kurve aufgetragen ist. 


M- KREIS 


N- KREIS 
N¢1 


Fig. 2 


Nyquist-Diagramm mit M-N-Kreisen. 


Diese Darstellungsart wird oft als Nyquist-Diagramm bezeichnet. Wenn die 
Gangkurve des aufgeschnittenen Regelkreises bekannt ist, so kann (9,/8.) (iw) — 
des geschlossenen Regelkreises mit Gleichung (9) bestimmt werden. Gleichung 


(9) kann auch graphisch gelést werden, wie in Figur 2 angedeutet ist. 
Es ist 


K Gy4|=OB, (18) 
1+ KG,4|=CB, (19) 


da C die Koordinaten (—1, 0) hat. 
Damit wird das Amplitudenverhialtnis 


(20) 


3 a) |= 
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Die Phasenverschiebung < (%,/9,) (i w) ist bestimmt durch: 
Oy F 
B=9—y= 52 (io), (21) 


wie aus Figur 2 ersichtlich. Es ist nicht schwierig, zu zeigen, dass die Kurven 

fiir verschiedene Werte von 
i 
Oo 


é 


(tw) | = const = M (22) 
und 
tg x $8 (iw) = const = N (23) 


im Polardiagramm von Figur 2 Kreise sind [4]. Ein M- und N-Kreis ist dort 
eingezeichnet. Ohne auf die Ableitung der Formeln fiir die M-N-Kreise hier 
naher einzugehen, seien der Vollstandigkeit halber die Gleichungen zur Be- 
stimmung der Radien und Zentren der Kreise in Tabelle 1 zusammengestellt. 


Tabelle 1 
| M-Kreise | N-Kreise 
? MM? 1 il 
Kreiszentrum | |- Wen’ 0| | [- 2° sn | 
Kreisradius een it : VN? a! 
| oe = — N 


Falls in einem Polardiagramm wie in Figur 2 viele M-N-Kreise eingetragen 
werden, so kann direkt der Frequenzgang (i,/,) (i w) des geschlossenen Regel- 
kreises bei gegebener Gangkurve K G;;,,) abgelesen werden. In Figur 2 tangiert 
K Gi») den M-Kreis im Punkte D. Eine Tangente zwischen der Gangkurve 
K Gi) und einem M-Kreise bedeutet ein Maximum oder Minimum im Ampli- 
tudenverhiltnis | (#,/0,) (¢ @)| des geschlossenen Regelkreises. Bei gegebener 
Gangkurve K G;;,) ist sofort ersichtlich, ob es sich um ein Maximum oder 
Minimum handelt. In Figur 2 zum Beispiel ist Punkt D ein Maximum. 


4. Stabilitat eines Regelsystems 


Eine Stabilitatsuntersuchung, sei es analytisch, graphisch oder durch 
Versuch, gehért zu den ersten Untersuchungen, die an einem Regelsystem 
unternommen werden. Es ist daher von Nutzen, ein Kriterium zur Hand zu 
haben, welches iiber die Stabilitat einer Regelung ohne grosse Mithe Auskunft 
gibt. Speziell gut geeignet bei Verwendung der Gangkurve K G;;,,) ist das 
Stabilitatskriterium von Nyquist. Die Ableitung des Kriteriums, entwickelt 
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fiir gegengekoppelte Verstarker, soll hier nicht weiter verfolgt werden, da dies 
den Rahmen dieses Aufsatzes iiberschreiten wiirde [5], [6], [7]. Die Formulierung 
des Nyquist-Kriteriums, angewendet auf einen einfachen Regelkreis fiir die_ 
Gangkurve K G;;,.), lautet: 

Wenn die Gangkurve —K G;;~) den Punkt (+1, 0) umschliesst, ist das 
System unstabil, andernfalls ist es stabil. 

Hierbei durchlauft die Frequenz m die Werte von —oo bis +00. Wir 
brauchen allerdings.nur die Gangkurve fiir positive Frequenzen aufzutragen, 


KURVE 1: STABIL wr Be Sa 
me 
KURVE TI: UNSTABIL w 


Fig. 3 
Stabilitatskriterium von NyguistT. 


da K G,_;,) durch Spiegelung von K G,, ;,,) an der reellen Achse erhalten wird. 
In Figur 3 stellt Kurve J eine stabile und Kurve JJ eine unstabile Regelung dar. 

Wenn die Gangkurve —K G;;,.. den Punkt (+1, 0) umschliesst, dann um- 
schliesst die Gangkurve K G;;,,. den Punkt (—1, 0). Das heisst, das Nyquist- 
Kriterium kann entweder auf den Punkt (+1, 0) angewendet werden, falls die 
Gangkurve —K G;;,,, aufgezeichnet wird, oder aber auf den Punkt (—1, 0) 
wenn man die Gangkurve +K G;;,) verwendet. Da in unsern Gleichungen 
immer + G;;,, auftritt, ist es unnétig, —K Gi», Zu bestimmen, und das 
Stabilitatskritertum wird deshalb in diesem Aufsatz durchwegs auf den Punkt 
(—1, 0) angewendet. 


—— 


5. Qualitative Abschitzung der Antwortkurve des geschlossenen 
Regelkreises mittels der Gangkurve K Gi.) 


Die Antwortkurve ist definiert als die zeitliche Anderung des Ausganges 
bei einer plétzlichen Anderung des Einganges eines Regelelements. In diesem 
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Abschnitt soll besonders die Antwortkurve des geschlossenen Regelkreises 
untersucht werden. In Figur 4 sind verschiedene Antwortkurven aufgezeichnet, 
wie sie haufig in der Praxis angetroffen werden. 

Bei der Synthese von Regelsystemen wird der Antwortkurve eine bestimmte 
Form vorgeschrieben. Die Antwortkurve ist durch Gleichung (12) mathema- 
tisch beschrieben. Der Zusammenhang zwischen der Anderung eines Regelele- 
mentes und der Verbesserung einer.bestimmten Wurzel der charakteristischen 
Gleichung (13) ist mit Ausnahme einfacher Regelkreise sehr untibersichtlich. 
Es ist daher schwierig, mittels der Antwortkurve rechnerisch oder graphisch 


EINGANG oh 


es 2 


* ii 


AuscaNG wh 


KURVE I : SCHWACHE DAMPFUNG 
KURVE I? KRITISCHEDAMPFUNG 


KURVE OI: STARKE DOAMPFUNG 


= 


0 t 
Fig. 4 
Antwortkurven fiir ein System zweiter Ordnung. 


eine optimale Regeleinstellung zu finden. Andererseits leistet die Antwortkurve 
sehr wertvolle Dienste bei der Optimaleinstellung von Regelungen im Betrieb 
oder bei der Bestimmung von Gangkurven aus Versuchsresultaten. Ist die 
charakteristische Gleichung von héherer Ordnung als 4, was in praktischen 
Beispielen fast immer der Fall ist, so wird der Arbeitsaufwand zur Bestimmung — 
der Wurzeln und deren Ubertragung auf den Entwurf von Regelelementen 
sehr gross. 

Es ist deshalb wiinschenswert, Wege zu finden, wo mittels der Gangkurve 
K Gj) direkt das Verhalten der Antwortkurve des geschlossenen Regelkreises 
bestimmt werden kann. Durch das Kriterium der M-N-Kreise sind wir sofort 
in der Lage, die Gangkurve (04/0) ( w) des geschlossenen Regelkreises, bei 
gegebener Gangkurve K Gj; ..), 2U bestimmen. In Figur 5 sind fiir den geschlos- 
senen Regelkreis, Amplitudenverhaltnis | (P4/P.) (@ w)| und Phasenverschie- 
bung @, separat aufgetragen. 

Eine Spitze im Amplitudenverhaltnis | (9./9, (¢ w) | deutet im allgemeinen 
auf eine komplexe Wurzel der charakteristischen Gleichung hin, wobei der 
imaginare Teil der Wurzel ungefahr der Frequenz entspricht, an welcher die 


ZAMP III/20 
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Spitze auftritt. Die Hohe der Spitze ist bestimmt durch den Realteil der kom- 
plexen Wurzel. Sollen nun die reellen und komplexen Wurzeln der charakte- 
ristischen Gleichung grosse Realteile haben, dann miissen die Spitzen von 
| (#,/3,) (¢ w) | begrenzt werden und sollten bei hohen Frequenzen auftreten. 
Hoch ist selbstverstandlich ein relativer Begriff und bedeutet hier: hoch ver- 
glichen mit dem Frequenzbereich, fiir welchen die Regelung entworfen ist. 
Es hat sich nun in der Praxis erwiesen, dass bei Beschrankung des maximalen 
Amplitudenverhiltnisses M = | (#,/8,) (¢ w)| auf 1,3 die positiven Realteile der 
Wurzeln in der Regel geniigend — 
gross sind, um die verschiedenen — 
Eigenfrequenzen zufriedenstel-_ 
lend zu dimpfen. Dieses Krite- 
rium ist selbstverstandlich nur 
eine Faustregel. Doch leistet es 
sehr gute Dienste bei der Lésung 
von Regelproblemen wie den hier 
behandelten. Die Wahl des maxi- 
malen Amplitudenverhiltnisses M 
wird je nach dem Zweck, den die 
Regelung erfiillen soll, verschieden 
sein. Zum Beispiel wird M haufig 
gleich 1 gewahlt, wodurch die 
Regelung kritisch gedimpft wird. 
Da es haufig nicht méglich ist, die elektrischen, mechanischen oder hydrau- 
lischen Komponenten zum voraus genau zu bestimmen, sondern Unsicherheits- 
faktoren von 20 bis 30% gebrauchlich sind, ist eine genaue Bestimmung der 
einzelnen Regelparameter fiir optimales Arbeiten der Regelung nicht gerecht- 
fertigt. Damit erhalt das M-Kriterium erhéhte Bedeutung, denn es erlaubt im 
allgemeinen, die optimale Einstellung innerhalb 30% zu bestimmen. Somit 
sind wir jetzt in der Lage, direkt vom Frequenzgang des aufgeschnittenen Re- 
gelkreises, mittels der M-N-Kreise, auf das Verhalten der Antwortkurve des 
geschlossenen Regelkreises zu schliessen. Da Anderungen des Frequenzganges 
K G;;,.) leicht auf physikalische Parameter tibertragen werden kénnen, so ist 
die Wirkung einer Parameteranderung auf die Antwortkurve sofort erkennbar. 


Fig. 5 


Gangkurve eines geschlossenen Regelkreises. 


6. Die Bestimmung des Frequenzganges K GE w) 


Beim Entwurf einer Regelung sind immer bestimmte Elemente als gegeben 
zu betrachten. Es ist dann Aufgabe der Synthese, zusiitzliche Regelelemente 
so zu bestimmen, dass zusammen mit den Gegebenen eine vorgeschriebene 
Gangkurve K G;;,) erhalten wird. In einem gegebenen Regelsystem ist phy- 
sikalisch am einfachsten, den Verstarkungsfaktor K der Gangkurve K G;;,.) zu 
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andern. Bei einer elektrischen Regelung zum Beispiel kann durch Anderung 
eines Spannungsteilers K leicht variiert werden, wahrend bei hydraulischen, 
pneumatischen oder mechanischen Regelsystemen die Anderung von K mei- 
stens mit Andern von Schlitzen, Ubersetzungsverhaltnissen und Federn ver- 
bunden ist. Die linke Abbildung in Figur 6 zeigt die Wirkung verschiedener 
auf die Gangkurve K G in einem Polardiagramm. Wie zu ersehen ist, ergibt 
eine Anderung von K nur eine Verlangerung oder Verkiirzung des Vektors, 
ohne dessen Phasenlage p zu andern. 


Es soll jetzt K Gnicht imlinearen, son- AMPLITUDE |KGI 
dern im logarithmischen Massstabe, mit KG= i z nea 
Amplitudenverhaltnis und Phasenlage $ 
getrennt, aufgetragen werden. Es ist: , imo: 
RS Aa 
RO teks cic.) ey eo as 
g| ee) | oe g | (i 5) ae T I:K=l PHASENVERSCHIEBUNG gunci 
Wr sike al) 
und 
gpk Gio) sa We Gio) . (25) Fig. 6 
’ : ( Effekt einer Verstarkungsfaktoranderung K 
Da die Elektroingenieure diese Be- auf die Gangkurve KG(iw). 


ziehungen seit langem benutzen, werden 
wir auch die dort gebrauchliche Einheit, das Dezibel (db), verwenden. Damit 
wird das Amplitudenverhaltnis A in Dezibel: 


A = 201g |K Gyq)| = 201g K + 20 lg | Gia] - (26) 


Die Phasenlage von wird durch Gleichung (25) bestimmt [8]. Durch diese 
Transformation von einem linearen zu einem logarithmischen Massstab werden 
Multiplikationen und Divisionen zu Additionen und Subtraktionen vereinfacht. 
Eine Anderung des Verstarkungsfaktors K bedeutet jetzt lediglich eine Parallel- 
verschiebung von | K Gijo)| nach oben oder unten, wie aus dem rechten Teil 
von Figur 6 ersichtlich ist. 

Wie in Figur 2 gezeigt wurde, kénnen die Werte fiir konstantes Amplituden- 
verhaltnis und Phasenlage des geschlossenen Regelkreises im Polardiagramm 
als Kreise dargestellt werden. 

Fs lasst sich nun ein neues Diagramm bestimmen, fiir 


Ge fe 
M = const = a Fo), und N = tg X | @@) = const, 


c 
in welchem die Ordinate in Dezibel und die Abszisse als Phasenverschiebung 
in Grad aufgetragen ist. Dieses Diagramm wird als Nichols-Diagramm bezeich- 
net [5]. Wie aus Figur 12 ersichtlich, ist die Gangkurve K Gj) im Nichols- 


Diagramm bestimmt durch | K Gia) | in Dezibel, am rechten Rand aufgetragen, 


308 RupoiF GLaus ZAMP *_ 


und PRG jem) 2M obern Rand ablesbar. Um (#,/0,) (¢ m) zu finden, haben wir 


nur die Schnittpunkte von K G;;,,, mit den stark ausgezogenen Kurven heraus- 
zulesen. Wir erhalten damit das Amplitudenverhaltnis 


$* (i 0) ak aeKe in Dezibel 
und die Phasenverschiebung 


e [5° @@)|=9 [ase] in Grad. 


Die Bestimmung des grésstméglichen Verstarkungsfaktors K erfolgt so, dass 
die Gangkurve G;;,,) vertikal verschoben wird, bis sie eine vorgeschriebene 
M-Kurve beriihrt. AK ist dann die Ordinatendifferenz zwischen K Gio) und 
Gio). Das Stabilitaétskriterium von Nyguist auf das Nichols-Diagramm ange- 
wendet, erfordert fiir einen einfachen, stabilen Regelkreis, dass die Gangkurve 
unterhalb dem Punkte (0 db, —180°) vorbeigehe. 
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Beziehungeén zwischen den Stérungskoeffizienten 
der ballistischen St6rungsrechnung 


Von Kart EcKEL, Madrid 


In der ballistischen Stérungstheorie werden die Anderungen einer Flugbahn 
berechnet, die durch kleine Anderungen der Flugbahnparameter entstehen. Als 
solche hat STANKE!) insbesondere die Anfangsgeschwindigkeit v), den Erhohungs- 
winkel #), den Bodenwert T, der absoluten Temperatur und den ballistischen 
Beiwert c in Betracht gezogen. Bezeichnet?) man die horizontale Bahnkoordinate 
mit x, die vertikale mit z, die Flugzeit mit ¢ und fiihrt die Gréssen uw und 0 durch 
die Gleichungen 


= Ino, = nts (4+) 


ein, so kann man die Anderungen Au, Az, Ax, At der Flugbahnvariablen, die 
infolge von gleichzeitig auftretenden Stérungen 40), Au, Ac und AT, verur- 
sacht sind, aus folgenden linearen Gleichungen berechnen: 


Ac AT, 
Au = O49 AOy + yy, AM + By 5 + yy ce 
0 
Ac ANGE, 
Az = 9 A@y + H%q1 Aug + Be Rene A, a Pa 
0 
(1) 
Ac Ade 
AX = gq Ay + Og, AUy + Be + Ys coat 
0 
Ac AT, 
At = O49 AQ + oq, AMq + Piers ainaied z 
0 


Vergleichspunkte auf der gestorten und ungestérten Flugbahn sind dabei die 
Punkte mit gleichem 9. Die Stérungskoeffizienten «;;, B,; und y; hangen von O 
ab und werden durch Integration von Systemen linearer Differentialgleichungen 
erster Ordnung ermittelt, die man als Stérungsgleichungen bezeichnet. Jede 
Integration liefert die vier Stérungskoeffizienten einer Spalte der Gleichungen (1). 
Mit Hilfe der von STANKE gefundenen Beziehungen zwischen den Stérungskoeffi- 
zienten «;,, 6; und y;, die als erweiterte Darrieussche Formeln bezeichnet werden, 
kann eine dieser Integrationen erspart werden, wenn die entsprechenden beiden 
anderen durchgefiihrt sind. Im folgenden soll eine weitere einfache Beziehung 
zwischen den Stérungskoeffizienten abgeleitet werden, in denen auch die Koeffi- 
zienten «; ) enthalten sind. 


1) P. STANKE, Die Storung einer Flugbahn durch dussere Einflisse, Webrtechn. Mh. 42, 560 


(1938); 43, 35, 63 (1939). 
2) Die Bezeichnungen sind iibereinstimmend mit R. SANGER, Ballistische Stoérungstheorte 


(Birkhauser, Basel 1950) gewahlt. 
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Zu einem Punkt P der ungestérten Flugbahn mit den Bahnkoordinaten u,, 2, 
Xp, O, und ¢, geh6re die Temperatur 7, und die Luftdichte 9,. Berechnet man 


mit u,, O, als neuen Anfangswerten, mit 7, als Bodenwert und mit c’= ¢ @,/Q9 : 


: 
: 


als ballistischem Beiwert eine neue Flugbahn vom Koordinatenanfangspunkt O — 


aus, so stimmt diese mit dem Teil der alten Flugbahn tiberein, der auf P folgt. 
Es findet nur eine Parallelverschiebung dieses Bahnteiles statt, der den Punkt P 
in O iiberfiihrt. Insbesondere bleiben die Werte 9 in den entsprechenden Punkten 
bei der Parallelverschiebung erhalten. Deshalb kann man, wenn die Differenzen 


AQ,= O,— O, Auy=uUuz—UMm, AC=c’—c, AT=—T,—-T 


klein genug sind, die Abweichungen der parallel verschobenen Flugbahn mit 
Hilfe der St6rungstheorie berechnen, indem man diese Differenzen in das Glei- 
chungssystem (1) einsetzt. Die linken Seiten dieses Gleichungssystems sind aber 
schon bekannt. Aw ist in allen Bahnpunkten gleich Null, da sich uw bei der Parallel- 
verschiebung nicht andert, 4z und Ax sind die Komponenten des Verschiebungs- 
vektors PO, und Af? ist in allen Bahnpunkten gleich dem negativen Wert der 
Flugzeit ¢,. 

Man berechnet die Stérungen 40), Au , Ac und AT, durch Taylor-Entwick- 
lungen um den Punkt O nach Potenzen von f,. Es gilt fiir den Punkt P, wenn 
die Entwicklung nach dem Glied erster Ordnung abgebrochen und die Erd- 
beschleunigung g konstant angenommen wird 


a) oe (2) 


Entsprechend ergibt sich die Stérung Aw, aus der Reihenentwicklung 


Te ke LY] aoe g [b ; 


wenn 6 die Verzégerung durch den Luftwiderstand bedeutet. Die Anderung des 
ballistischen Beiwertes ist durch die Anderung der Luftdichte AQ = Op — Oo 
bedingt, und es gilt 


Ac _ Ag 
¢ 20 
Fiir das lineare Temperaturgesetz 
T=T%—Az 
erhalt man 
Ac 1 (42) dz 
os sae (ST AE oa 
c Cn \ zi aw at.) 1 
Sod Of RAAT URED ee 
ae ( Fe :) Vp Sin By t, + ++ (4) 
1 SAS: 
= —Uy ao Sin Oy ft +s, 


"Ra 


worin F& die Gaskonstante fiir Luft bedeutet. 


a et i oe a ptm ae 


Vol. III, 1952 Kurze Mitteilungen — Brief Reports - Communications breves Silt 


Fiir das gleiche Temperaturgesetz wird 
AT, 1 ( aI dz 
emai aiers%> 
Wt [A EN ie 


d 1 v : ane 
ae ( Bi 2] My Sin Og by 2° (5) 


= = Sin ge pei s = 
Setzt man die Stérungen (2) bis (5) in (1) ein, dividiert beide Seiten durch ¢ 
und geht zur Grenze ¢, > 0 iiber, so verschwinden die vernachlassigten Glieder 
hoéherer Ordnung. Der Verschiebungsvektor PO geht in die Tangentenrichtung 
der ungestérten Flugbahn in O iiber und ist der Anfangsgeschwindigkeit v9 ent- 
gegengesetzt gerichtet. Daher wird 


: Az : 
lim ——=—v,sind?, lm mess —U, Cost . 
typ? 90  “p tp—0 “p 
Ferner gilt Az/t, = —1, so dass man insgesamt folgende Beziehungen fiir die 


Stoérungskoeffizienten erhalt: 
0 = Ag Mo + Ay % + B By ae Cy, 
Vp Sin By = Ag Xo + A,%, + Bp,+ C ye, | 
Vg COS Hy = Ag to + Ay %i + B Ag+ Cys, | 
Le Age + Alen + B pet Ove; 


worin die vier Konstanten durch 


4 o> Pos 4 24 tg bo in 9 
x Vo” eros ( ee zt a), 
(Bes ey, — = i Cele i 
i (5 2) sin By , . Vo Sin By 


gegeben sind. 

Die Richtigkeit dieser Formeln kann man durch Einsetzen in die St6rungs- 
gleichungen verifizieren. Man hatte auch umgekehrt an Stelle des hier gegebenen 
anschaulichen Beweises durch Einsetzen von (6) in die Stdrungsgleichungen mit 
anschliessendem Koeffizientenvergleich die Koeffizienten Ay, Ay, B, C und die 
linken Seiten von (6) berechnen k6énnen. 

Zur Bestimmung aller Stérungskoeffizienten wird man praktisch zunachst 
durch Integration des einfachsten homogenen Systems der Stérungsgleichungen 
die Koeffizienten «;, und durch Integration des einfachsten inhomogenen Systems 
die Koeffizienten B; berechnen. Mit Hilfe der erweiterten Darrieusschen Formeln 
erhalt man dann die Koeffizienten y; und durch obige Beziehungen (6) die Koeffi- 
zienten aj. 


Résumé 


Les déviations des variables de la trajectoire u, 7, z et? peuvent étre calculées 
suivant STANKE a partir des équations (1), comme fonction des perturbations 
AO, Atty, Ac et AT, lorsque Von connait les coefficients de perturbation corres- 
pondants a, %, B; et 7: Introduisant dans ces équations comme perturbations 
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les déviations qu’un point P de la trajectoire normale avec le temps de vol 7#, 
possede en comparaison avec le point d’origine, les cétés gauches de (1) sont 
égaux a 0, —z,, —%,, —t,. En développant les perturbations autour du point © 
d’origine en puissances de ¢,, en divisant par ¢, et en se rapprochant de la limite { 
t, — 9, on obtient les relations linéaires (6) pour les coefficients de perturbation. 
Us forment un complément aux formules de DARRIEUs et ses élargissements de 
STANKE, ‘ 


(Eingegangen: 8. April 1952.) 


Automatische Rechenplanfertigung bei programmégesteuerten 
Rechenmaschinen?) 


Von HEINZ RUTISHAUSER, Ziirich?) 


Sn a I el ow 


Es ist bekannt, dass die Aufstellung des Rechenprogramms (Befehlsreihe) fiir 
ein bestimmtes Problem oft nicht nur eine zeitraubende Arbeit, sondern auch eine 
Quelle von sehr stérenden Fehlern sein kann. Im Gegensatz zu andern Vorschla- 
gen zur Vereinfachung der Planfertigung [insbesondere K. ZusE%)], welche die 
Konstruktion spezieller Gerate in Betracht ziehen, zeigt der Verfasser einen Weg, 
den Rechenplan fiir eine bestimmte Formel oder Formelgruppe durch die Rechen- 
maschine selbst «ausrechnen» zu lassen, das heisst die Rechenmaschine als ihr _ 
eigenes Planfertigungsgerat zu beniitzen. 

Zu diesem Zweck werden die Formeln zunachst arithmetisiert, das heisst auf — 
numerische Form gebracht, indem man jedem «Element» (das sind Klammer- ‘ 
zeichen, Operationszeichen und Operanden) eine Zahl zuordnet. Diese Zahlen 
sind fiir Operations- und Klammerzeichen ein fiir allemal festgelegt (zum Bei- 
spiel 10000 fiir eine 6ffnmende Klammer, 60000 fiir ein Multiplikationszeichen | 


le oe 


usw.). Fiir einen Operanden dagegen ist es die Nummer der Zelle, die man zur 
Speicherung desselben vorgesehen hat. 

Die so erhaltene Zahlfolge charakterisiert eine Formel oder Formelgruppe — 
volistandig, so dass es méglich sein muss, daraus den Rechenplan zu berechnen. 
Es wird eine Methode angegeben, wie man diese Zahlfolge in endlich vielen ; 
Schritten abbaut und gleichzeitig die gesuchte Befehlsreihe aufbaut. : 

Interessant wird die Sache allerdings erst bei zyklischen Problemen, also — 
bei Formeln, in denen laufende Indizes vorkommen. Hier hat man librigens : 
zwei Moglichkeiten: man kann namlich entweder einen zyklischen Plan herstellen — 
(der mit einem Minimum an Befehlen auskommt) oder einen sogenannten «ge- | 
streckten» Plan, welcher jede auszufiihrende Operation ohne Riicksicht auf die 
zyklische Natur der Aufgabe genau einmal enthalt. Gestreckte Rechenplane sind 

1) Autoreferat tiber einen vom Verfasser an der GaMM.-Tagung 1951 in Freiburg i. Br. gehal- 
tenen Vortrag, dessen Inhalt demnichst als Mitteilung aus dem Institut fiir angewandte Mathe- 
matik an der ETH. veréffentlicht werden wird. ; 

Uber Begriffe und Bezeichnungen des programmgesteuerten Rechnens vergleiche: H. Ruris- 
HAUSER, A. SPEISER, E, STIEFEL, Programmegesteuerte digitale Rechengerdte, Mitteilung Nr. 2 aus 
dem Institut fiir angewandte Mathematik der ETH. (Birkhauser, Basel 1951); ZAMP 1, 297 und 
339 (1950); 2, 1 und 63 (1951). 

2) Institut fiir angewandte Mathematik der ETH. 


8) K. Zusr, Uber den allgemeinen Plankalkiil als M ittel zur Formulierung schematisch-kombina- 
tiver Aufgaben, Arch. Math. 1, 441-449 (1948) 
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zwar sehr lang, rechnen aber dasselbe Problem wesentlich schneller durch, so 
dass man mit Vorteil davon Gebrauch machen wird, wenn man einmal der 
manuellen Herstellung enthoben ist. 

Der Verfasser gibt ein Rezept fiir die automatische Herstellung des zyklischen 
oder gestreckten Rechenplans fiir jedes Problem, das sich in ahnlicher Weise 
wie das folgende Beispiel der Matrizeninversion formulieren lasst: 

Sei a>, die zu invertierende Matrix, dann lautet die Inversionsvorschrift, 
welche tiber die Zwischenmatrizen a, (j = 1, 2,...,2%—1) die gesuchte Inverse 


als a’, liefert?) 


Lee) = all) aes 
Eis 2 


lehoapy;>——= ae ee i 


> ne ar ; 
Fir k= 2(1)n:0,, Xa, = a@),,_4- 
EM — 02) (1 70): 
Fir & = 1: Soo ag = at 


ner oe Sat a 
2Furlk = 2 (1) 224.5 = (ai, x al, rt = yy 4: 


Eine spezielle Anwendung: 
Es seien Differentialgleichungen zweiter Ordnung, 


y” = f(#, VV’), 


numerisch zu integrieren. Normalerweise macht man das so, dass man einen 
permanenten Plan fiir «Integration einer gewohnlichen Diffefentialgleichung 
zweiter Ordnung» hat und daran einen Unterplan fiir die spezielle Formel 
f(x, v, y’) anfiigt; den letzteren kann man aber nach den beschriebenen Methoden 
berechnen, man muss nur /(#, y, y’) in numerische Form kleiden, wozu das Plan- 
fertigungsgerat dienen soll. Wenn damit auch kein grosser Zeitgewinn verbunden 
sein mag, so eliminiert man doch durch das Ausschalten der manuellen Rechen- 
planfertigung eine Fehlerquelle. 


Summary 


It is well known that planning and coding of problems for an electronic com- 
puter is very often a laborious task and a source of errors. To simplify the coding 
(programming), various devices have been constructed, but the author believes 
that the device most suited to assist the coding personnel is the computing 
machine itself. Methods to «compute» the program for a given set of formulae 
have been worked out by the author and are described in detail in a booklet to be 
published very soon (see footnote 1 on page 312). 


(Eingegangen: 1. Mai 1952.) 


1) Diese Vorschrift basiert auf der sogenannten Jordanschen Elimination, doch werden dort 
die Zwischenmatrizen in der Regel nicht berechnet. Man beachte ferner, dass in der folgenden Vor- 
schrift die zu berechnenden Gréssen immer rechts vom Gleichheitszeichen auftreten. 

2) 4 =a(b) c ist die gebrauchliche Abkiirzung fiir: «i durchlauft die Werte a, a + b,a+ 2b, 


. usw., bis c». 
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Uber die numerische Berechnung der Eigenwerte von 
Sakulargleichungen 


Von HEINZ VOETTER, Wien?) 


Eine n-reihige Sakulardeterminante (1)*) ist eigentlich 


a,— ¥ Ay Un | 
91 Ag — *¥ Gal) Seg (1) 
a3 A39 Azy 


eine algebraische Gleichung m-ten Grades, die in entwickelter Form (2) bekannt- 
lich lautet: 
n 
1 
i 


ie 1)i a zara = 0. (2) 


A) bedeutet dabei einen ¢-zeiligen Hauptminor der Determinante (1). Aus (2) 
k6nnen mit Hilfe des Horner-Schemas in Verbindung mit der Regula Falsi die 
Eigenwerte der Gleichung (1) rasch berechnet werden. Schon bei Sakularglei- 
chungen fiinften Grades wird aber die Ausfiihrung der Formel (2) sehr langwierig, 
da die Anzahl der zu berechnenden Determinanten stark anwichst. 

Es ist unseres Wissens noch nirgends darauf hingewiesen, dass die Verwandlung 
der Determinante (1) in die algebraische Form (2) durch «Triangulierung» mit 
verhaltnismassig geringem Aufwand méglich ist. Das Verfahren sei an einer 
Sakulardeterminante vierten Grades gezeigt. Es lasst sich ohne weiteres auf 
Determinanten beliebigen Grades erweitern. 


Oy eo as AN3 AN4 | 41 
Ao, Agg — ¥ Aggy Ao4 Ag2 
D=| ay A39 433 — ¥  Agg A43 (3) 
M1 M1 + Ay1 MQ + Ag1 Ag 


Aq Ago A43 Reg Fo aN A 
Aq) 


Wir wollen durch Kombination eine neue Zeile der Determinante (3) erhalten, 
bei der das Glied in der ersten Spalte eine Null wird und die in allen anderen 
Spalten mit Ausnahme der letzten keine Glieder mit der Unbekannten x enthilt. 
Die einzelnen Elemente in den Zeilen in (3) werden mit den rechts stehenden 
Zahlen multipliziert und dann addiert. Die Multiplikatoren der Zeilen, die (mit 
Ausnahme der letzten Spalte) ein x enthalten, sind aus der untersten Zeile zu 
entnehmen. Die unterste Zeile wird jeweils mit der Unbekannten multipliziert, 
wodurch letztere in allen Spalten der neuen Zeile mit Ausnahme der letzten 
wegfallt. Der konstante Summand des Multiplikators der letzten Zeile ist so 
bestimmt, dass das Glied in der ersten Spalte der neuen Zeile wegfallt. Wir er- 
halten damit eine neue Determinante (4). Gegeniiber (3) ist die oberste Zeile weg- 
gelassen, und als letzte die auf eben beschriebene Weise erhaltene angeschrieben. 


1) Institut fiir Verfahrenstechnik, Technische Hochschule Wien. 
*) In der mathematischen Literatur werden nur die symmetrischen (4;; = a;;) Determinanten 
vom Typ (1) als Sakulargleichungen bezeichnet und die ubrigen als «charakteristische Gleichungen». 


—- 


a ee 
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Ay, Ag9 — ¥ Ags Ao4 by 

31 Age a33 — © ~ Aza by 
oe Qo, 01 + Ae, 
= 21 y 

D Ay 42 43 Aggy = 22 1 31 <2 A (4) 
gy 
Gani a a 
Ome: b, PRA by ie Woy, = eT ieee a A 
1 


Das Verfahren wird sinngemiass fortgesetzt. Die nachste neue Zeile soll in den 
ersten beiden Spalten den Wert Null besitzen und in der dritten kein Glied mit 
« haben, was durch die angeschriebenen Multiplikatoren erreicht wird. A besei- 
tigt das Glied der ersten Spalte und der konstante Multiplikator der vierten Zeile 
das Glied der zweiten Spalte. Nach nochmaliger Wiederholung des Vorganges 
erhalt man die Determinante (5). 


Chey Urey ies CEs a2 

OQ 0) lay “Ona Ga 

OM OMMINC IMIG pi CAs CAE at = 

OD © O Chae Gh oP ee Uh eee Chee > Be 


(5) 


Das rechts unten stehende Polynom in (5) reprasentiert die entwickelte alge- 
braische Form der Sakulargleichung. Bei der Durchfiihrung dieses Verfahrens 
braucht nicht wie’in (3) und (4) usw. nach jedem Schritt die vollstandige Deter- 
minante angeschrieben zu werden. Die neuen Zeilen werden der Reihe nach unter 
die Ausgangsdeterminante gesetzt. Als Beispiel fiir eine rationelle Schreibweise 
sei eine Sakulargleichung vierten Grades ausgefiuhrt. 


2,858 — x\— 0,377 0,152 

2,846 — ¥—0,157 0,856 2,4514 | 

S561 Oa 0787 0,288 0,1644 0,5775 

0,152 0,856 0,288 0,650 — mv + 0,3495 | 5,6273 3,1116 
2,4514 0,1644 0.3857 + 0,3005 —x2|x —4,7582 | —0,9011 

0,5775 2,1246 —6,6714¥-+5,0587 x2—x8 |v —2,0492 

2 5303 + 12,4133 v—16,1366 x2 + 7,1079 43— x4 

Brahe ayaa? 10sL ee 


Als Probe fiir die richtige Durchfiihrung dieser Operation kann der Koeffizient 
der zweithéchsten Potenz der Unbekannten dienen. Dieser muss nach (2) gleich 
sein der Summe der Elemente der Glieder der Hauptdiagonale von (1). ; 

Die Ermittlung der Multiplikatoren und alle Zwischenrechnungen sind mit 
einer normalen Rechenmaschine ohne weitere Zwischennotizen mdoglich, was bei 
SAkulardeterminanten auch héheren Grades eine rasche Durchfiihrung der Rech- 


nung gewahrleistet. 


Summary 


The calculation of the eigenvalues of secular-equations is a difficult and 
tedious task. Even with a five-ordered determinant it is very time-consuming 
to evaluate the minors in the usual way. The purpose of this paper is to describe 
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a simple method to transform the secular determinant into its algebraic form. 
By consequent combination of rows of the determinant all elements which are 
standing under the main-diagonal are eliminated, so that the value of the de- 
terminant is represented only by the product of the diagonal elements. This 
method can be used with the normal calculating-machine. The tabular scheme 
of the calculation also makes it possible for a quick execution of secular-equa- 
tions of higher degrees. 

(Eingegangen: 7. Februar 1952.) 


Varia — Miscellaneous — Divers 


Organisation fiir angewandte Physik und Mathematik in der Schweiz 


Schon vor einiger Zeit ist der Wunsch aufgetaucht, es méchten die in unserem 
Lande auf dem Gebiet der angewandten Physik und Mathematik Tatigen sich in 


ee nm 


a 


einer passenden Organisation zusammenfinden. In der Tat ist das Bediirfnis — 


nach einer ungezwungenen Darlegung und Diskussion der erhaltenen Forschungs- 
ergebnisse recht gross. Wahrend auf dem Gebiete der reinen Physik, der Chemie, 
der reinen Mathematik usw. entsprechende Vereinigungen bestehen, sind die 


es 


«Angewandten» nirgends zu Hause; ihre Themata kénnen sich in den Rahmen : 


der dort gehaltenen Vortrage nicht recht einfiigen. Angesichts des sehr betracht- 
lichen Umfanges der angewandten Forschung ist das ein wenig befriedigender 
Zustand. — Eine naheliegende Lésung ware die Schaffung einer eigenen Gesell- 
schaft fiir angewandte Physik und Mathematik gewesen. Allein von Seiten der 
reinen Physik wurden Bedenken gedussert, ob eine Zersplitterung der Physik 
(bzw. Mathematik) in einen reinen und einen angewandten Teil nicht doch grosse 
Nachteile brachte. Der aus idealen und praktischen Griinden erwiinschte Kontakt 


wiirde noch geringer werden als bisher. Da diesem Argument ein grOsseres Gewicht | 


beizumessen ist, wurde von einigen interessierten Herren in Besprechungen eine 
Vereinbarung mit der Schweizer Physikalischen Gesellschaft (SPG.) vorgeschlagen 
— in dem Sinne, dass die SPG. an ihren beiden Tagungen im Friihjahr und im 
Herbst (letztere zusammen mit der Naturforschenden Gesellschaft) in ihren Tages- 
programmen den angewandten Richtungen angemessen Rechnung tragt und eben- 
so bei der Zusammensetzung des Vorstandes. Darnach wire die Griindung einer 
besonderen Gesellschaft unnétig; an Stelle dessen wiirde die einfache Mitglied- 
schaft zur SPG. geniigen. 

Um gewissermassen experimentell zu erproben, wie eine solche Zusammen- 
arbeit funktionieren kann, wurde anlasslich der Friihjahrstagung der SPG., am 


3. Mai 1952, in Baden, parallel mit den kernphysikalischen und allgemein-physi- — 


kalischen Vortragen eine Serie von 13 Kurzreferaten iiber angewandte Themata 
angesetzt. Dank dem Entgegenkommen der Firma Brown, Boveri AG. und einer 
sehr zahlreichen Beteiligung erwies sich dieser Versuch als in erfreulichem Masse 
gelungen. Kurze Inhaltsangaben der Referate werden in der ZAMP erscheinen, 
einzelne Referate auch in extenso. 

An der geschaftlichen Sitzung der SPG. wurde der Anregung zugestimmt, 
und damit ist die Méglichkeit von regelmassigen Tagungen geschaffen. — Als 
zweiter Sekretar und als Vertreter der angewandten Richtung wurde Dr. P. pr 
HALLER, Gebr. Sulzer AG., Winterthur, gewahlt. 

Die nachste Versammlung der SPG. findet vom 23. bis 25. August 1952 in 
Bern statt. J. ACKERET 
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Theoretische Elektrotechnik. IIL: Grundziige dey Theorie elektrischer 
Maschinen. Von K. KUHLMANN (Birkhauser, Basel 1951). 547 S., 328 Abb.; 
sFr. 74.90. 


Das vorliegende Buch von Professor K. KUHLMANN, der dritte Band seines 
Sammelwerkes tiber die theoretische Elektrotechnik, in dem in grossen Ziigen 
seine an der Eidgendssischen Technischen Hochschule gehaltenen Vorlesungen 
iiber dieses Gebiet niedergelegt sind, behandelt einige der wichtigsten Probleme 
der elektrischen Maschine. Im Gegensatz zu vielen ahnlichen Biichern werden 
nicht die einzelnen Maschinentypen in sinngemasser Reihenfolge besprochen; 
hier wird das einzelne Problem in seiner Vollstandigkeit durchgearbeitet und auf 
seine Anwendung bei den verschiedenen Maschinen hingewiesen. Die Behandlung 
der Probleme wird vorzugsweise in analytischer Form durchgefiihrt, was natur- 
-gemass manchmal zu etwas langlichen Formeln und zu Gréssensymbolen mit 
einigen Indizes fiihrt. Der klare und logische Aufbau der Uberlegungen und die 
saubere Darstellung der Entwicklung lassen den Leser doch verhaltnismassig 
leicht folgen. 

Der Stoff des Buches gliedert sich in acht Kapitel, deren Inhalt im folgenden 
kurz angedeutet sei. Nach einem kurzen, einige allgemeine Festlegungen umfas- 
senden Kapitel A folgen im Kapitel B unter der Bezeichnung der Wirkungsweise 
der elektrischen Maschinen die grundlegenden Feststellungen tiber Fluss und 
EMK. Kapitel C befasst sich mit den Wicklungen, fiir deren Schaltgesetze ana- 
lytische Ausdriicke angegeben werden. Bei dem auf Seite 65 beschriebenen Bet- 
‘spiel einer doppelten Schleifenwicklung hatte vielleicht darauf hingewiesen wer- 
den kénnen, dass durch andere Wahl des Nutenschrittes die beiden Spannungs- 
polygone zur Deckung gebracht werden und dadurch Ausgleichstr6me vermieden 
werden kénnen. Das folgende Kapitel D enthalt die umfassende analytische Be- 
handlung der Erzeugung magnetischer Felder verschiedenster Formen und Wick- 
lungsausfiihrung. Kapitel E ist der Induktivitat gewidmet. Hier werden die wich- 
tigen Rechenformeln fiir Eigen-, gegenseitige und Streuinduktivitat abgeleitet. 
Die auf der Induktivitat beruhende Berechnung der induzierten EMK. bildet den 
Inhalt des folgenden Kapitels F, wobei auch der Einfluss der héhern Harmoni- 
schen in seiner Vollstandigkeit behandelt ist. Die magnetische Energie und das 
Drehmoment finden ihre Durcharbeitung in Kapitel G. Das folgende Kapitel H 
ist als Anhang bezeichnet, indem hier einige Sonderprobleme Aufnahme gefunden 
haben, die zum Teil nur lose mit der Maschine zusammenhangen. Der Leser findet 
hier die Durcharbeitung des Kommutationsproblems als induktiver Vorgang. 
Im weitern hat hier die analytische Behandlung des Problems der Stromverdran- 
gung in runden und rechteckigen Leitern und im Eisen Aufnahme gefunden, 
ebenso die Berechnung der Eisenverluste. Ein letzter Abschnitt in diesem grossen 
Kapitel orientiert iiber das Rechnen mit symmetrischen Komponenten und dessen 
Anwendung auf die Vorgange in den verschiedenen Leitersystemen. 

Es darf ruhig behauptet werden, dass Professor KUHLMANN mit diesem Buch 
ein erstklassiges Werk geschrieben hat, das in sehr klarer Form unter Bentitzung 
einer konsequenten Schreibweise die wichtigsten elektromagnetischen Zusam- 
menhange der elektrischen Maschine festlegt. Die zahlreichen, den Text erlau- 
ternden Zeichnungen sind ausserst sorgfaltig ausgefiihrt und sehr instruktiv. Auch 
die Ausstattung, welche der Verlag Birkhauser dem Buche zuteil werden liess, 
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verdient volle Anerkennung und hilft mit, dem Buche ein Niveau zu geben, das — 
nur selten von ahnlichen Biichern erreicht wurde. Das Werk sei einem recht 

zahlreichen Leserkreis als etwas vom besten, das in diesem Gebiet bisher geboten ~ 
wurde, zum Studium warmstens empfohlen. E. Diinner 


my So 


Ecoulements des fluides compressibles. Par R. SAvER (Béranger, Paris 
et Liége 1951). 307 p., 150 fig.; fr.fr. 3900.—. 


Le professeur SAUER est connu d’un public particuliérement vaste du fait 
qu’il écrivit en 1943 un remarquable petit livre sur la théorie des écoulements 
compressibles stationnaires (Einfiihrung in die theorvetische Gasdynamik). Cet 
ouvrage, remanié et enrichi de la matiére de tout un cours sur les phénoménes 
non stationnaires, a été publié en francais: La traduction du manuscrit allemand 
est due a un aérodynamicien autorisé, M. CARRIERE. 

Apres un élégant exposé des équations générales pour les écoulements sans 
friction ni conductivité, des solutions simples sont obtenues par linéarvisation des 
dites équations. Il faut remarquer que ces solutions, qui sont approximatives pour 
un gaz, sont exactes dans le cas du mouvement non stationnaive d'un liquide dans 
un tube, et trouvent ainsi les applications les plus variées. Etudiant en théorie 
linéaire la répartition de pression sur un corps de vévolution dans un écoulement 1 
stationnaire coaxial, l’auteur n’arrive a déceler aucune influence de la compressi- 
bilité. Nous pensons que ce n’est pas l’approximation linéaire de la compressibilité 
qui est insuffisante, mais plut6t la maniéve dont l’écoulement incompressible de— 
base est lui-méme linéarisé!). 

A partir de l’approximation linéaire, et par des démonstrations plutét géo- 
métriques, le lecteur arrive sans peine aux solutions rigoureuses dans le domaine — 
supersonique. La méthode des caractéristiques, dont il s’agit ici, s’applique aussi — 
aux écoulements non stationnaires et aux écoulements non isentropiques. Dans 
le premier cas (unidimensionnel) une interprétation graphique trés simple, et dans 
le second cas des conseils pratiques pour le calcul numérique, dénotent le souci 
de rendre la théorie immédiatement utilisable. C’est grace a cet esprit pratique 
que l’auteur, tout mathématicien qu'il est, sait mettre en évidence les méthodes 
les mieux adaptées aux problémes de la technique. Aussi pour I’étude des flux 
stationnaires subsoniques ou transsoniques, les procédés lents mais stirs des dé-_ 
veloppements en série et des approximations successives sont-ils exposés avec le © 
plus de détail. 

Dans la derniére partie du livre, consacrée aux ondes de choc et aux écoule- 
ments — stationnaires ou non — qui en contiennent, les lois régissant les ondes 
de détonation intéresseront plus spécialement le thermodynamicien. Un tableau 
des grandeurs avant et aprés le choc nous semble unique dans la «littérature 
numérique» actuelle du fait que la grandeur d’entrée est le taux de compression 
au travers du choc, paramétre commode puisque invariant a l’égard des trans- 
lations de lobservateur. 

D’un paragraphe a l’autre et d’un bout a l’autre du livre s’étend un véritable 
réseau de comparaisons, soulignant les similitudes ou les différences entre les 
nombreux cas étudiés. Les résumés placés en téte des chapitres, ainsi que les 
phrases d’orientation judicieusement insérées dans le texte, sont les meilleurs 
guides au travers de cet ouvrage si riche. B. Chaix 


a 


——e 


1) Icile passage a la limite d'un rayon nul est tel qu’il ne subsiste point de longueur a laquelle 
on puisse rapporter ce rayon avant d’en prendre le logarithme. 
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Die zweidimensionale Laplace-Transformation. Von D. VoELKER und 
G. Dortscu (Verlag Birkhauser, Basel 1950). 259 S., Fr. 39.-. 


Das vorliegende Werk ist ein Lehrbuch der Transformation 


CoO Cc 


& (F(x, y)} = fu, v) = f fe-™-99 F(x, y) dx dy, 
0 0 


welche bei Funktionen von zwei Variablen dasselbe leistet wie die gewéhnliche 
Laplace-Transformation bei Funktionen einer Variablen und deshalb fiir die 
Lésung von partiellen Differentialgleichungen eine Rolle spielt. Das Buch — zur 
Zeit das einzige iiber dieses Gebiet — widmet sich deshalb auch weniger der Theo- 
rie, die im wesentlichen der klassischen Theorie analog ist, sondern mehr den An- 
wendungen auf partielle Differentialgleichungen, deren Grundgebiet eine Viertel- 
ebene ist. 

Dabei wird nun sehr sch6n gezeigt, wie man solche Gleichungen — etwa die 
Warmeleitungsgleichung und die Wellengleichung — lost und vor allem wie sich 
dariiber hinaus die Vertraglichkeitsbedingungen mit Hilfe der zweidimensionalen 
Laplace-Transformation von selbst ergeben. Es werden auch Anwendungen auf 
Systeme von partiellen Differentialgleichungen sowie auf partielle Differential- 
gleichungen mit drei unabhangigen Variablen behandelt. Die zweite Halfte des 
Buches enthalt auf zirka 100 Seiten ausftihrliche Tabellen von Korrespondenzen 
fiir die zweidimensionale Laplace-Transformation /(u, v)<—>F(#, v), die man fiir 
die Anwendungen sehr schatzen wird. AL, Rutishauser. 


Numerische Behandlung von Differentialgleichungen. Von L. CoLLAtz 
(Springer-Verlag, Berlin 1951). 458 S., 110 Abb.;.DM 45.—. 


Das'‘vorliegende Werk gibt eine vollstandige Ubersicht tiber die numerischen 
Methoden zur Lésung von Differentialgleichungen. Es behandelt sowohl Anfangs- 
als auch Randwertprobleme bei gew6hnlichen wie bei partiellen Differential- 
gleichungen und wird jedem, der sich in der Praxis mit solchen Aufgaben befassen 
muss, eine wertvolle Hilfe sein, nicht zuletzt auch wegen seines ausgezeichneten 
systematischen Aufbaus. 

Der erste Abschnitt behandelt Anfangswertprobleme bei gew6hnlichen, Dif- 
ferentialgleichungen, wobei die Methoden vom Typus RuNGE-Kutta und 
Apams-BAsHFOoRT im Mittelpunkt stehen, wahrend der zweite Abschnitt vor 
allem die Anwendung der Differenzenverfahren und der Verfahren vom Ritzschen 
Typus auf Rand- und Eigenwertprobleme zeigt. Der dritte Abschnitt beschreibt 
die Lésung parabolischer und hyperbolischer partieller Differentialgleichungen 
vor allem mit Differenzenverfahren, wobei insbesondere auch die Fragen der 
Fehlerfortpflanzung und Stabilitat der Verfahren eingehend besprochen werden. 
Schliesslich wird in einem vierten Abschnitt gezeigt, wie man elliptische Diffe- 
rentialgleichungen lést. Insbesondere wird auch hier auf die Relaxationsmethode, 
das Ritzsche Verfahren und andere Minimalprinzipien eingegangen. Ein ftinfter 
Abschnitt ist den Funktional- und Integralgleichungen gewidmet. Den Schluss 
bildet eine Sammlung von Formeln zur Integration von Differentialgleichungen, 
von Ausdriicken fiir das Differenzenverfahren, von Eulerschen Ausdriicken und 
vor allem von Operatoren fiir die Relaxationsrechnung. 

Dabei ist der Text laufend durch numerische Beispiele erganzt, die zeigen, 
dass der Verfasser eine grosse praktische Erfahrung besitzt. Auch ist hervorzu- 
heben, dass iiberall grosse Sorgfalt auf Fehlerabschatzungen verwendet wurde. 

AH, Rutishauser. 
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Integraltafeln, Sammlung unbestimmter Integrale elementarer Funk- 
tionen. Von W. MEYER zUR CAPELLEN (Springer-Verlag, Berlin 1950). VIII + 
292 S.; DM. 36.—. ‘ 

Das Buch stellt eine umfangreiche Sammlung von Integralen dar.- Diese sind 
nach Integranden geordnet: Integrale rationaler Funktionen (48 Seiten), Inte- 
grale algebraischer Funktionen (50 Seiten), ausserdem elliptische Integrale (38 Sei- 
ten). Von besonderer praktischer Wichtigkeit ist der grosse Abschnitt tiber trans- 
zendente, insbesondere trigonometrische und hyperbolische Integrale (80 Seiten). 
Ihm folgt eine Auswahl von Integralen, deren Integranden sowohl algebraische 
als auch transzendente Bestandteile enthalten. 

Fiir zahlreiche Typen von haufig gebrauchten Integralen werden allgemeine 
Gesetze oder Rekursionsformeln angegeben, mit deren Hilfe die Zuriickfiihrung 
auf explizite angefiihrte Integrale méglich ist. Dadurch wird die Verwendbarkeit 
der Tafeln natiirlich betrachtlich erh6éht. 

Diese auch drucktechnisch sehr gut ausgeriistete Integraltafel wird zweifellos 
dem Praktiker ausgezeichnete Dienste leisten; doch wird auch der reine Mathe- 
matiker gerne zu diesem Hilfsmittel greifen. H. Rutishauser 


Einfiihrung in die Funktionentheorie. Von L. BIEBERBACH (Verlag fiir 
Wissenschaft und Fachbuch, Bielefeld 1951). 220 S., 43 Abb.; DM. 12.60. 


Cet ouvrage est une excellente introduction aux parties de la théorie des 
fonctions analytiques d’uhe variable qui doivent étre connues de tout étudiant 
en mathématiques et en physique (il serait, a ce point de vue, désirable qu’une 
prochaine édition contienne un bref exposé de l’application de la théorie aux 
théoremes d’existence des équations différentielles, spécialement des équations 
différentielles linéaires), Grace a la limpidité de son style, il peut étre reeommandé 
aussi a lingénieur qui veut utiliser les méthodes de la théorie des fonctions en 
théorie du potentiel et en hydrodynamique. L’auteur ne suppose du lecteur que 
la connaissance des éléments du calcul différentiel et intégral. 

M. Plancherel 


Mechanische Schwingungen. Von J. P. DEN HARtoG und G. MESMER 
(Springer-Verlag, Berlin 1952), 427 S., 299 Abb.; DM 42.-. 


Die englische Sprache hat sich im Westen heute so weit durchgesetzt, dass 
-man sich fragen kann, ob sich eine Ubertragung englischer Lehrbiicher in die 
deutsche Sprache noch lohnt. Wenn ja, dann sicher im vorliegenden Falle, denn 
das 1947 in 3, Auflage erschienene Werk von J. P. DEN Hartoc, Mechanical 
Vibrations, ist einzig in seiner lebendigen und anschaulichen Darstellungsart 
sowie in seiner im Einzelnen knappen, im Ganzen umfassenden Schau. Nachdem 
schon die erste Auflage (1934), aus einem Kurs des Verfassers an der Design 
School der Westinghouse Company hervorgegangen, von G. MESMER ins Deutsche 
tibertragen worden ist, war eine deutsche Bearbeitung der stark erweiterten 
dritten Auflage (1947) fallig, und wenn sich wiederum G. MEesMER dieser Aufgabe 
unterzogen hat, so hat er sich um so mehr Dank und Anerkennung verdient, als 
er nicht eine blosse Ubersetzung, sondern eine Uberarbeitung des englischen 
Textes gibt. Der Originaltext erscheint darin in vielen Einzelheiten verbessert 
und erganzt, und insbesondere das Literaturverzeichnis hat durch Erweiterung 
unter vermehrter Beriicksichtigung europdischer Autoren erheblich an Wert 
gewonnen., H. Ziegler. 
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